User's  Manual  for  the 
Scanning  Fast  Field  Program 
(SCAFFIP) 

General  Version  1.0 


by  John  M.  Noble 
Dave  Marlin 

Battlefield  Environment  Directorate 


19950217  001 


Approved  for  public  release ;  distribution  is  unlimited 


NOTICES 


Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of 
the  Army  position,  unless  so  designated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in  this  report  is  not  to 
be  construed  as  official  Government  indorsement  or  approval  of  commercial 
products  or  services  referenced  herein. 


Destruction  Notice 

When  this  document  is  no  longer  needed,  destroy  it  by  any  method  that  will 
prevent  disclosure  of  its  contents  or  reconstruction  of  the  document. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


3.  REPORT  TYPE  AND  DATES  COVERED 


1.  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 

User’s  Manual  for  the  Scanning  Fast  Field  Program  (SCAFFIP) 

General  Version  1.0  _ _ 


6.  AUTHOR(S) 

John  M.  Noble,  Dave  Marlin 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


U.S.  Army  Research  Laboratory 
Battlefield  Environment  Directorate 
ATTN:  AMSRL-BE-S 


9.  SPONSORING /MONITORING  AGEUCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-545 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


ARL-TR-545 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

The  Scanning  Fast  Field  Program  (SCAFFIP)  is  an  atmospheric  acoustic  propagation  model 
incorporating  many  of  the  effects  on  the  environment  on  the  sound  field  such  as  geometrical 
spreading,  refraction,  diffraction,  molecular  absorption,  and  complex  ground  impedance.  SCAFFIP 
provides  attenuation  levels  with  range  and  azimuth  or  sound  pressure  levels  in  dB  (re.  20  pPa)  with 
range  and  azimuth  for  a  given  geometry,  frequency,  and  meteorological  profile.  The  meteorological 
profile  and  the  geometry  provides  the  model  the  ability  to  calculate  the  sound  speed  profile.  The 
geometry  profile  is  required  because  the  angular  dependence  of  the  sound  speed  on  the  wind 
direction  is  relative  to  the  direction  of  propagation.  This  model  works  over  a  flat-earth  and  a  non- 
turbulent  atmosphere.  Even  with  these  restrictions,  the  model  performs  very  well  for  many 
scenarios.  The  model  contains  a  friendly  user  interface  requiring  a  minimum  amount  of  information 
to  run  the  model.  There  are  also  flags  that  can  be  set  to  obtain  more  detailed  information. 


:  14.  SUBJECT  TERMS 


ii  FFP,  acoustic,  model,  atmospheric,  Fast-Field  Program 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17  SECURITY  CLASSIFICATION  i  18.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  ,  OF  THIS  PAGE  OF  ABSTRACT 


Unclassified 


NSN  7540-01-280-5500 


j  SAR  _ j 

Standard  Form  2S5  (Rev  2-89; 

Prescribe  by  AN$i  StcJ  Z39'18  • 

298-102  1 


Contents 


1.  Introduction 


2.  Theory  .  5 

2.1  Speed  of  Sound  .  5 

2.2  FFP  .  8 

2.3  Absorption  of  Sound  in  the  Air .  16 

2.4  Complex  Ground  Impedance  .  19 

3.  Operations .  23 

3. 1  Purpose  of  User  Interface .  23 

3.2  Format  of  Input  Data  Files  .  23 

3.2.1  Debug  .  23 

3.2.2  Weather .  24 

3.2.3  Geometry  .  25 

3.2.4  Frequency .  26 

3.2.5  Ground  .  26 

3.3  Explanation  of  User  Interface .  27 

3.4  How  to  Run  SCAFFIP  .  27 

3.5  Format  of  the  Output .  28 

References .  29 


Acronyms  and  Abbreviations 


31 


Appendices 

Appendix  A.  Source  Code  for  SCAFFIP  General  Version  1.0 .  33 

Appendix  B.  Example  Cases  .  83 

Appendix  C.  Debug  Features .  109 

Appendix  D.  User  Developed  Interface  Requirements  for  SCAFFIP .  113 


Distribution 


.  *099.8940^  .Fpr.  ....... 

HTIS  GRA&I  Ef* 

DUG  TAB  □ 

Unannounced  □ 

Tua*  1 nnt.1  fln 

Dtat^ibut  ton  JMfc, 

Availability  Gpdea 

r— - 

Met. 

V' 

Avail  and/or 
Special 

117 


Figures 


1 .  Diagram  of  geometry  definition  . 

2.  Layering  of  the  atmosphere  by  the  FFP . 

3.  Transmission  line  analog  used  by  the  FFP . 

4.  Log-log  plot  of  sound-absorption  coefficient  versus  frequency 

for  sound  in  air  at  20  °C  at  1  atm  pressure  and  with  a  water-vapor 
fraction  h  of  4.676  x  10"3  (Rh  =  20  %)  . 

Tables 


1.  Transmission  line  analog  relationships  to  acoustic  fields  . 

2.  Estimation  of  0  and  a  for  various  ground  surfaces  given 

n'  =  0.750  and  Sf  =  0.875  . 


2 


1.  Introduction 


The  Scanning  Fast  Field  Program  (SCAFFIP)  is  based  on  the  Fast  Field 
Program  (FFP)  with  the  added  ability  to  scan  multiple  azimuths  to  predict  the 
propagation  conditions  about  the  location  of  a  sensor.  SCAFFIP  makes  a 
prediction  of  the  acoustic  propagation  conditions  based  on  spherical  spreading, 
molecular  absorption,  refraction,  acoustically  complex  ground  impedance,  and 
diffraction  over  benign  terrain. 

The  FFP  is  a  one-way  solution  to  the  acoustic-wave  equation  originally 
developed  for  underwater  sound  propagation  predictions.  [1,2]  Raspet  et 
al.  [3]  and  Lee  et  al.  [4]  adapted  the  FFP  to  propagation  in  the  atmosphere. 
The  FFP  developed  by  Raspet  et  al.  used  a  propagation  matrix  formulation. 
If  each  layer  in  the  atmosphere  is  viewed  as  an  optical  device,  a  matrix  for 
each  layer  in  the  atmosphere  can  be  constructed.  Multiplying  each  matrix 
together  results  in  a  new  matrix  representing  how  an  acoustic  signal  will  be 
affected  as  it  propagates  through  the  atmosphere.  Next,  a  Bessel  Function 
Transform  of  the  problem  is  taken  with  respect  to  range.  After  the  solution  is 
calculated,  an  inverse  transform  is  performed  to  arrive  at  the  final  solution. 

The  software  package  that  comes  with  SCAFFIP  is  an  integrated  set  of 
algorithms  for  running  the  acoustic  propagation  model  on  any  computer 
platform  containing  an  ANSI  version  of  a  C  compiler.  There  is  an  option  of 
using  the  nongraphical  user  interface  that  comes  with  the  package  for  extracting 
the  propagation  model  and  incorporating  it  into  the  user  interface  environment 
(see  appendix  D).  Example  test  cases  (in  appendix  B)  show  the  structure  of  the 
input  data  files  used  by  SCAFFIP  and  give  a  basis  to  work  with  to  investigate 
the  installation  of  the  files.  The  software  package  is  available  on  a  variety  of 
computer  media.  A  list  of  possible  media  is  4-mm  Digital  Audio  Tape  (TAR 
or  CPIO  format),  0.25-in.  cartridge  tape  (TAR  or  CPIO  format),  8-mm 
Exabyte  tape  (TAR  or  CPIO  format),  MS-DOS  3.5-in.  high-  or  low-density 
diskettes,  or  0.5-in.  9-track  magnetic  tape  (TAR  or  CPIO  format). 
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2.  Theory 


2.1  Speed  of  Sound 


Meteorological  phenomena  can  have  a  significant  effect  on  the  received  sound 
field.  Some  of  the  meteorological  variables  affecting  the  speed  of  sound  in  air 
are:  pressure,  temperature,  vector  wind  speed,  and  humidity.  To  observe  the 
effect  of  each  meteorological  variable,  consider  each  one  independently  and 
examine  the  equation  for  the  speed  of  sound  in  air.  The  value  of  c,  according 
to  Laplace’s  adiabatic  assumption  for  air  as  an  ideal  gas,  [5]  is 


c  (T)  = 


N 


Y  RT 
M 


(1) 


where 

7  =  the  ratio  of  specific  heats 

R  =  the  universal  gas  constant  equal  to  8314.16  J/(kg  K) 

M  =  the  molecular  weight. 

The  presence  of  water  molecules  alters  the  sound  speed  by  lowering  7  and 
decreasing  M.  The  decrease  in  M  dominates  so  that  the  overall  effect  of 
increasing  humidity  is  an  increasing  sound  speed.  These  changes  can  be 
quantified  as 


Y  = 


7  +h 
5  +h 


(2) 


and 


M  =  29-11  h 


(3) 


where 

h  =  the  fraction  of  water  molecules  in  air. 
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As  the  amount  of  water  in  the  atmosphere  increases,  the  molecular  weight  of 
a  parcel  of  air  will  decrease  because  the  molecular  weight  of  a  water  molecule 
is  less  than  diatomic  nitrogen.  This  molecular  weight  effect  will  try  to  increase 
the  sound  speed  as  the  fraction  of  water  molecules  in  the  air  increases.  To 
calculate  the  fraction  of  water  molecules  in  air,  the  Goff-Gratch  equation 
equation  (4),  must  be  used  to  first  calculate  the  partial  pressure  of  saturate 
water  vapor,  Psat  at  temperature  T. 


(pJ\ 

=  10.79586 

1  - 

/  T  \ 

-  5.02808  log  10 

1  T ' 

T 

kJ 

■ 

{  T 

roij 

+  0.42873 x  10-3(l0«-sll-‘r”'m)  -  l)  -25195983 

where 

T01  =  273.16  K 

po  =  1  atm  or  the  reference  pressure. 

After  the  value  for  Psat  is  determined,  the  fraction  of  water  molecules  in  air  can 
be  calculated  using  the  following  relationship: 

,  10  -2(RH)Psat  (5) 


where 

RH  =  the  relative  humidity  in  percent 
P  =  is  the  pressure  in  atmospheres. 

The  magnitude  of  the  dependence  of  the  sound  speed  on  humidity  is  not 
obvious.  To  understand  the  degree  of  the  effect  of  humidity  on  sound  speed, 
consider  a  particular  case.  At  20  °C,  the  difference  in  sound  speed  between 
0  and  100  percent  humidity  is  2  m/s.  A  fluctuation  in  the  humidity  of  this 
amount  is  very  unlikely.  If  the  variation  in  humidity  is  reduced  to  a  change  of 
50  to  100  percent,  the  change  in  the  sound  speed  is  only  1  m/s.  Therefore,  the 


variation  of  sound  speed  caused  by  changes  in  humidity  should  always  be  much 
less  than  1  m/s.  Generally,  humidity  fluctuations  can  be  ignored. 

The  effect  of  the  wind  speed  on  the  speed  of  sound  is  a  vector  relation.  The 
effective  sound  speed  is  calculated  using 

Cgff  =  c(T)  +  u  •  cos  (0w-rc  -6*)  (6) 

where 

c(T)  =  the  speed  of  sound  in  the  absence  of  wind  at  temperature  T 
u  =  the  magnitude  of  the  horizontal  wind  speed 

0R  =  the  bearing  of  the  receiver  from  the  source 

0W  =  the  direction  from  which  the  wind  blows 
0w  -  7i"  =  the  direction  the  wind  is  blowing  (figure  1). 

All  directions  are  relative  to  north. 

The  sound  speed  will  also  vary  with  height  because  the  sound  speed  is  a 
function  of  temperature  and  vector  wind  speed.  This  variation  will  cause  the 
acoustic  wave  to  be  refracted  as  it  propagates  through  the  atmosphere.  The 
degree  of  refraction  the  acoustic  wave  undergoes  is  related  to  the  sound  speed 
gradient  present  in  the  atmosphere.  If  the  sound  speed  increases  with  height, 
the  acoustic  wave  will  be  refracted  downwards.  If  the  sound  speed  decreases 
with  height,  the  acoustic  wave  will  be  refracted  upwards. 
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2.2  FFP 

The  propagation  of  sound  from  a  point  source  located  at  the  origin  is  given  by 
the  classical  wave  equation 


V2p  -  —  =  -4u6(x,y,z) 

c2  dt2 


(7) 


where  8  represents  a  delta  function  source  of  unit  strength.  For  simple 
harmonic  motion,  equation  (7)  becomes  the  Helmholtz  equation 


V2p  +  k2p  =  -4%b(x,y,z) 


(8) 
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where 


k  =  is  the  wavenumber,  co/c,  in  which  c  —  the  sound  speed,  co  —  the 
angular  frequency. 

For  the  FFP,  k  and  c  are  restricted  to  vary  only  in  the  z-direction. 

Transforming  equation  (8)  into  cylindrical  coordinates  and  assuming  azimuthal 
symmetry,  the  Helmholtz  equation  becomes 


J*E  +  i  Op  +  +  e-p  =  ~-6(r)Hz-z,) 

dr2  r  dr  3z\p  dz)  r 


(9) 


where  the  source  is  located  at  r  =  0  and  z  —  zs  and  p  is  the  density  of  the 
medium. 


The  atmosphere  is  viewed  as  a  series  of  constant  sound  speed  layers  for  the 
FFP  (as  shown  in  figure  2).  The  layers  in  the  atmosphere  are  bounded  on  top 
and  bottom  by  complex  impedance  surfaces.  The  top  boundary  is  typically 
modeled  as  an  infinite  half-space  with  constant  parameters.  At  the  bottom 
boundary,  the  atmospheric  layer  adjoins  a  partially  absorbing  surface  that  can 
be  represented  by  the  complex  acoustical  impedance  of  the  ground. 


To  reduce  the  dimensionality  of  equation  (9),  a  zero-order  Hankel  transform 
is  applied  with  respect  to  the  rangej/ariable  r.  This  gives  the  transform  pair: 


p(K,z)  =  f  p(r,z)J0(Kr)rdr 

O 


(10) 


p(r,z)  =  f  p(K,z)J0(Kr)KdK. 

O 
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Applying  the  first  transform  to  equation  (9)  results  in 


+  [k\z)  -  k 2]p  =  -26(z-zs). 
dz 2 


LayerN  ImpedanceZ* 


Figure  2.  Layering  of  the  atmosphere  by  the  FFP. 


10 


This  equation  can  be  decomposed  into 


v  i  dp 
u  = - — 

*'  «  Po  dz 


(12) 


— —  ^-2-  =  -  — —  [  k2(z)  -  k2  ]  p  +  — —  5(z-zs)  (13) 

dz2  <»P 0  i“Po 


mean  air  density 
uz  =  the  transformed  particle  velocity  in  the  z-direction. 

The  delta  function  is  the  source  term  producing  a  discontinuity  in  uz  at  a  height 
zs.  To  solve  the  resulting  equations,  Lee  et  al.  [4]  used  an  analogy  to  a 
transmission  line,  which  results  from  the  form  of  the  transformed  equation. 
The  transformed  equations  have  a  very  similar  form  to  the  telegrapher’s 
equations  of  electrical  transmission  line  theory.  From  the  analogous 
telegrapher’s  equations, 

—  =  - ZI(z )  (14) 

dz 


and 

duz 

dz 

where 
Po  =  the 


and 

dl(z)dz  =  -Y(z)V(z)  +  I0&(z-zs)  (15) 


there  is  a  similar  form  of  the  equations  if  the  shunt  admittance  Y(z)  is  made  a 
function  of  z  and  the  series  impedance  Z  is  made  a  constant  with  a  current 
source  at  z  =  zs.  The  equivalent  transmission  line  configuration  to  figure  2  is 
shown  in  figure  3. 
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Figure  3.  Transmission  line  analog  used  by  the  FFP. 

Using  the  Lee  et  al.  analogy,  the  acoustic  problem  can  be  arranged  so  that  a 
solution  can  be  calculated.  The  analogy  is  made  by  representmg  each  layer  m 
the  atmosphere  by  an  element  in  a  transmission  line  with  a  certain 
characteristic  admittance  and  attenuation  constant.  The  admittance  is  defined 
as  one  over  the  impedance  or  one  over  the  sum  of  the  resistance  plus  the 
reactance  of  the  electrical  element.  The  admittance  of  the  element  causes  the 
voltage  running  through  it  to  be  attenuated  and  a  shift  m  the  phase  of  the 
signal.  Viewing  the  analogy  from  the  perspective  of  the  acoustic  wave 
problem,  as  the  acoustic  wave  is  propagated  through  the  atmosphere,  it 
undergoes  losses  and  phase  shifting  caused  by  refraction  and  spreading  of  the 
acoustic  wave.  The  equivalence  between  the  transmission  line  and  the  acoustic 
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wave  problem  is  illustrated  in  table  1 .  This  equivalence  can  be  carried  out  for 
each  layer  of  the  atmosphere;  thus,  constructing  a  transmission  line.  The  top 
and  bottom  boundaries  in  the  atmosphere  become  loading  admittance  elements 
on  each  end  of  the  transmission  line.  The  problem  now  has  been  converted 
from  determining  the  voltage  in  the  transmission  line  to  a  point  in  the  line. 
This  is  a  well-known  process  in  electrical  engineering. 

Table  1.  Transmission  line  analog  relationships  to  acoustic  fields 


Acoustic  Wave  Transmission  Line 


Fields  Pressure  p(k,z) 

Velocity  6)  (k,z) 

Attenuation  Constant  Y  =  Vk2  -  (°Voc)2 
Characteristic  Admittance  Yc  =  -  y/cop 


Voltage  V(z) 
Current  I(z) 

y= VzTT 
yc=Vy7z 
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The  zero-order  Bessel  function  in  equation  (10)  can  be  expanded  using  Hankel 
functions: 


(16) 


The  Hankel  functions  can  be  represented  as  an  incoming  and  outgoing  acoustic 
wave.  The  FFP  is  designed  to  model  radially  outgoing  acoustic  waves.  This 
allows  for  the  first  Hankel  function  to  be  suppressed  with  the  additional 
argument  that  the  incoming  acoustic  wave  will  not  contribute  significantly  to 
the  final  result.  The  asymptotic  expansion  of  the  second  Hankel  function  for 
large  arguments  is  as  follows: 


aJV) 


e-i( Kr  -  it/4) 


xr  >  1 


(17) 


making  the  problem  easier  to  handle.  The  important  contributions  from  the 
integrand  of  the  inverse  transform,  equation  (10),  comes  from  the  area  where 
k  ~  k0.  Substituting  equation  (17)  into  the  inverse  transform  of  equation  (10) 
and  taking  the  far-field  approximation,  the  acoustic  pressure  equation  can  be 
written  as 


P(r,zr) 


^  +  f  p  (K,zr)  e'iKr  \/ic  dK. 
y/27i :r  o 


(18) 


To  perform  the  calculation  on  a  computer,  the  continuous  integral  must  be 
replaced  by  a  numerical  integral  over  discrete  values  of  k.  Applying  this  to 
equation  (18)  yields 


p(r,zT ) 


(1  +  0 


JV-l 


g-i  (2it nrn/N) 


(19) 


where 


Ak  = 

N-  1 
r 

m  =  — 

A  r 


K  =  bAk 

n 

a  2u 

A  r  =  - . 

AT  Ak 


(20) 


The  term  xmax  comes  from  the  property  of  the  integrand  of  equation  (10)  and 
only  has  significant  contributions  in  a  finite  range  of  k,  allowing  the  summation 
to  be  terminated  at  a  finite  number  of  terms. 

A  problem  in  the  derivation  is  one  of  the  complex  numerical  integrations.  The 
function  being  integrated  contains  branch  points  and  poles  on  the  real  axis. 
Because  of  the  nature  of  branch  points  and  poles,  the  integration  being 
performed  must  not  include  any  of  these  points  for  a  correct  solution.  To 
avoid  these  problems,  the  current  FFP  uses  what  is  called  extra  loss  in  the 
calculations.  The  mathematical  result  of  using  this  extra  loss  is  to  move  the 
numerical  integration  off  the  real  axis.  The  effect  of  the  extra  loss  is  removed 
from  the  solution  in  an  approximate  manner  by  multiplying  the  computed 
pressure  by  the  term  exp(ar),  where  a  is  the  extra  loss  attenuation  constant  in 
Np/m.  The  proper  choice  of  the  artificial  attenuation  is  essential  if  meaningful 
results  are  to  be  obtained  from  the  code. 

Another  problem  with  the  numerical  integration  is  the  number  of  points  N  used 
in  the  summation  of  equation  (19).  A  lower  bound  on  the  number  of  points 
required  in  the  summation  is 

N.  =  [6]  (21) 

min  _  L  J 


However  in  most  cases,  this  lower  bound  is  too  large  to  perform  one  single 
Fast  Fourier  Transform  (FFT).  The  problem  is  if  there  are  enough  points  to 
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sufficiently  sample  the  wavenumber  space  in  the  numerical  integration.  If 
equation  (19)  is  rewritten  in  the  form 

--iQnnmlN)  (22) 

Une 

n- 0 


N- 1 


this  summation  [6]  can  be  rewritten  so  that  the  single  summation  is  rearranged 
to  a  double  summation  of  the  form 


p 


b-  1 


N'-l  , 

„-i[2Tt(b-l)mlp]  V  a  p-i[2nnm/(pN )] 

e  Urn+(i>-l)iV/C 

n-0 


(23) 


where 

N '  =  N/p 

p  =  an  integer  larger  than  1. 

p  is  the  number  of  panels  that  the  original  FFT  has  been  divided  into.  Each 
panel  contains  N'  points.  This  technique  allows  the  FFT  to  be  performed  in 
order  to  calculate  the  acoustic  pressure  with  range  at  the  height  of  the  receiver. 

2.3  Absorption  of  Sound  in  the  Air 


Losses  in  the  medium  are  basically  caused  by  viscosity ,  heat  conduction,  and 
molecular  exchanges  of  energy.  In  the  nineteenth  century,  only  the 
mechanisms  of  viscosity  and  heat  conduction  were  suspected  of  causing 
dissipation  of  sound;  therefore,  they  are  presently  referred  to  as  classical 

absorption. 

In  classical  absorption,  if  the  effect  of  absorption  is  represented  by  a  factor  e 
where  r  is  the  distance  of  propagation,  then  the  attenuation  coefficient  acl 
caused  by  viscosity  and  heat  conduction  is  given  by  equation  (17)  from 
Physical  Acoustics  XVII  [7]: 
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(24) 


ucl  =  5.578  x  10  '9 


T!T0  f 
r +110.4  p/p0‘ 


The  units  of  acl  in  equation  (24)  are  nepers/meter,  where 

P0  =  the  reference  pressure  of  1.01325  X  105  Npm"2  (1  atm) 

P  =  pressure  in  Npm'2 
T0  =  the  reference  temperature  of  293.15  K 
T  =  temperature  in  K 
f  =  frequency  (Hz) 

In  molecular  absorption,  energy  exchanges  at  the  molecular  level  include 
rotational  and  vibrational  modes.  Analysis  of  the  rotational  mode  shows  that 
the  representative  attenuation  coefficient  is  proportional  to  o:cl,  the  classical 
attenuation  coefficient: 


^2  =  4.16  e‘16-8T'l/3 


a 


cl 


(25) 


when  293  K  <  T  <  690  K. 

For  frequencies  below  10  MHz,  it  has  been  demonstrated  that  energy  losses 
caused  by  classical  and  molecular  absorption  are  additive. 


A  simplified  empirical  form  of  the  equation  can  be  written  as 


a 


cr 


1.83  x  10 -11 


PIP0  ’ 


(27) 


which  is  correct  within  2  percent  for  213  K  <  T  <  373  K. 
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The  vibrational  mode  of  absorption  should  also  be  considered.  Because  the 
atmosphere  is  composed  mostly  of  nitrogen  and  oxygen,  each  will  contribute 
an  attenuation  coefficient,  where  j  stands  for  either  oxygen  or  nitrogen.  The 
symbols  are  defined  as  follows: 


cc 


vib,  j 


4  pXj 
35  c 


V2 

\T, 


(28) 


Xj  =  the  mole  fraction  of  air  component  considered  (0.20948  for  oxygen  and 
0.78084  for  nitrogen) 

q.  =  the  characteristic  vibrational  temperature  (2239.1  K  for  oxygen  and 
3352.0  K  for  nitrogen) 

c  =  the  speed  of  sound  at  temperature  T  (m/s). 


The  fr  j  are  the  individual  relaxation  frequencies  for  oxygen  and  nitrogen.  The 
computation  of  these  frequencies  depends  on  the  relative  humidity  and 
atmospheric  pressure.  The  relaxation  frequencies  are  given  by: 


P_ 

P0 


1 24  +  4.04  x  10 Ah 


0.02  +  h  \ 
0.391  +  h) 


f  =- 

Jr,N  p 
o 


X  T 


£  (  9  +  280  n/rr-q 

rr  ' 


(29) 


The  total  attenuation  coefficient  is  then  the  sum  of  acr  and  avib  j.  Figure  4  is 
a  log-log  plot  [6]  of  total  attenuation  coefficients  for  T  =  20  °C  and 
Rh  =  20  percent.  Figure  4  shows  the  attenuation  caused  by  classical 
absorption,  vibrational  relaxation  of  nitrogen  and  oxygen,  and  the  total 
attenuation  coefficient  caused  by  the  sum  of  the  three  attenuation  mechanisms. 

The  attenuation  coefficient  a  is  proportional  to  the  square  of  the  frequency. 
As  the  frequency  doubles,  the  attenuation  will  quadruple.  The  attenuation  of 
the  sound  wave  caused  by  molecular  and  vibrational  absorption  is  very 
important  for  frequencies  over  250  Hz.  For  frequencies  below  250  Hz,  this 
attenuation  does  not  contribute  much  to  the  total  attenuation  of  the  sound  wave. 


20  °C  at  1  atm  pressure  and  with  a  water  vapor 
fraction  h  of  4.676  x  10'3  (Rh  =  20  %). 


2.4  Complex  Ground  Impedance 


There  are  several  models  available  for  calculating  the  complex  ground 
impedance.  The  impedance  model  used  in  SCAFFIP  is  the  Four  Parameter 
Model  of  Attenborough.  [8]  In  this  impedance  model,  the  complex  normalized 
characteristic  impedance  of  the  ground  is  calculated  using 


4  q2  .Sfa 

— 2_  +  i— — 

3Q  co  pp 


(30) 
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where 


q2  =  fi-n' 

Sf  =  the  pore  shape  factor  ratio 

fl  =  the  porosity  of  the  ground 

a  =  the  flow  resistivity  of  the  ground  (mks)  rayls 

co  =  the  angular  frequency 

p0  =  the  density  of  air  (1.2  kg/m3) 

kb  =  the  normalized  wave  number. 


The  normalized  wave  number  kb  is  computed  from 


kb  “  \fyQ 


'A 

,3 


Y  ~  1 
Y 


+  i 


Sjo 


«P0 


1/2 


(31) 


where 

7  =  the  ratio  of  specific  heats  (equation  (2)) 

Npr  =  the  Prandtl  number  (0.724). 

The  parameters  Sf,  0,  p,  and  n'  are  normally  varied  until  agreement  is  reached 
between  impedance  measurements  and  the  impedance  model  is  achieved. 
However,  this  method  of  determining  the  four  parameters  cannot  always  be 
completed  if  time  or  resources  are  lacking.  Table  2  provides  rough  estimates 
of  the  values  for  the  four  parameters,  given  some  general  descriptions  of  a 
variety  of  ground  surfaces,  to  aid  the  user  when  the  four  parameters  are 
unknown. 

Some  ground  surfaces  have  a  layered  structure  resulting  from  the  gradual 
deposition  of  material  over  a  soil  base  (a  layer  of  snow  over  frozen  ground, 
decomposition  of  organic  material  over  clay  or  sandy  soil,  or  a  well  plowed 
pasture  over  clay  or  harder  packed  soil.)  An  effective  impedance  Z(d)  can  be 
calculated  for  a  semi-infinite  layer  of  impedance  Z2  covered  by  a  layer,  depth 
d,  of  another  material  of  impedance  Zt.  The  effective  impedance  is  given  by 


(32) 


Z(d)  = 


iZ1tan(kbd) 
iZ2tan(kbd)  1 


where 


kb  =  the  bulk  propagation  constant  in  the  top  layer  of  the  ground. 

The  parameters  Z,,  Z2,  and  kb  are  calculated  using  equations  (30)  and  (31) 
from  Attenborough’s  Four  Parameter  impedance  model. 
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Table  2.  Estimation  of  ft  and  a  for  various  ground  surfaces  given  n'  =  0.750 
and  Sf  =  0.875  _ 


Description  of  Surface 

ft 

alces) 

Dry  snow,  new  fallen  0.1  m  over 
about  0.4  m  old  snow 

0.850 

23 

Sugar  snow 

0.825 

48 

In  forest,  pine,  or  hemlock 

0.825 

48 

Grass:  rough  pasture,  airport,  public 
buildings,  etc. 

0.675 

330 

Roadside  dirt,  ill-defined,  small  rocks  up 
to  0.1  m  mesh 

0.575 

960 

Sandy  silt,  hard-packed  by  vehicles 

0.475 

3470 

"Clean"  limestone  chips,  thick  layer 
(0.01  to  0.025  m  mesh) 

0.425 

6470 

Old  dirt  roadway,  fine  stones  (0.5  m  mesh) 
interstices  filled 

0.400 

7500 

Earth,  exposed  and  rain-packed 

0.350 

17100 

Quarry  dust,  fine,  very  hard-packed 

0.300 

41700 

Asphalt,  sealed  by  dust  and  light  use 

0.250 

120000 
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3.  Operations 


3.1  Purpose  of  User  Interface 

The  user  interface  is  designed  to  run  the  model  with  numerous  scenarios  and 
a  minimum  of  modifications  to  the  input  files  between  each  run.  The  user 
interface  also  allows  the  model  to  be  run  interactively  or  in  a  batch  job  format. 
This  user  interface  will  be  standard  for  all  of  the  benign  terrain  acoustic 
propagation  models  released  by  the  U.S.  Army  Research  Laboratory  Battlefield 
Environment  Directorate  (ARL/BED).  This  allows  for  the  same  data  files  to 
be  used  for  future  benign  terrain  acoustic  propagation  models  obtained  from 
ARL/BED.  If  other  data  input  formats  are  required,  ARL/BED  can  assist  with 
the  modifications  to  the  current  user  interface.  Under  most  conditions, 
ARL/BED  cannot  guarantee  compatibility  with  future  upgrades  and  acoustic 
propagation  models  after  the  original  user  interface  has  been  modified. 

3.2  Format  of  Input  Data  Files 

There  are  five  input  files  required  by  SCAFFIP.  The  file  names  can  be 
anything  within  64  characters  in  length,  file  names  must  still  conform  to  the 
limitations  set  by  the  computer  system  being  used  because  SCAFFIP  asks  for 
the  name  of  each  file  as  it  runs.  Example  input  files  are  included  with  the 
discussion  of  test  cases  in  appendix  B.  The  five  categories  required  for  input 
are:  (1)  debug,  (2)  weather,  (3)  geometry,  (4)  frequency,  and  (5)  ground.  The 
syntax  of  each  of  the  five  files  is  given  in  this  section: 

3.2.1  Debug 

The  debug  file  contains  an  integer  code  that  allows  for  input  of  additional 
parameters  normally  fixed  within  the  program  and  output  of  additional  files. 


The  format  for  the  debug  file  follows: 


Code: 

0  Normal  operation 

1  Output  speed  of  sound  profile  (sound. o) 

2  Use  radians  for  wind  direction 
4  Source  centered  calculation 

8  Input  extra  loss  (extra. los) 

16  Input  number  of  points  and  panels  for  performing  FFT  (trans.frm) 

32  Output  panel  values  (npan.e) 

64  Output  ground  impedance  (imped. e) 

128  Output  wave  number  spectrum  (wavnum.e) 

256  Disable  smoothing 

The  value  for  code  can  be  one  of  the  numbers  listed  above  or  a  sum  of  any  two 
or  more  of  the  numbers  listed  above.  If  normal  operation  of  the  model  is 
desired,  use  a  zero  in  the  debug  file.  If  the  speed  of  sound  profile  and  wave 
number  spectrum  output  is  desired,  the  value  in  the  debug  file  would  be 
129  =  1  +  128.  Appendix  C  contains  the  file  formats  for  the  six  debug  files. 

3.2.2  Weather 

The  weather  file  contains  the  atmospheric  profile.  The  first  line  of  the  file  is 
an  integer.  The  rest  of  the  file  contains  floating  point  numbers. 

The  format  for  the  weather  file  follows: 


Nint 

Z(l) 

Z(2) 

T(l) 

T(2) 

Rh(l) 

Rh(2) 

P(l) 

P(2) 

Wvel(l) 

Wvel(2) 

Wdir(l) 

Wdir(2) 

Z(Nint) 

T(Nint) 

Rh(Nint) 

P(Nint) 

Wvel(Nint) 

Wdir(Nint) 
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where 


Nint  -  number  of  interfaces.  (There  is  no  limit;  however,  many 

interfaces  will  slow  the  program  down.  Try  and  limit  the 
number  of  interfaces  to  no  more  than  60.) 

Z(i)  -  height  of  the  ith  interface  in  meters,  also  Z(l)  =  0  and 

Z(i)  <  Z(i+1)  always. 

T(i)  -  temperature  in  Kelvin  at  height  Z(i). 

Rh(i)  -  relative  humidity  in  percent  at  height  Z(i). 

P(i)  -  pressure  in  atmosphere  at  height  Z(i). 

Wvel(i)  -  wind  speed  in  meters  per  second  at  height  Z(i). 

Wdir(i)  -  direction,  in  degrees,  from  which  the  wind  is  blowing 
relative  to  north  at  height  Z(i). 

3.2.3  Geometry 

The  geometry  file  contains  the  information  describing  the  geometry  of  problem. 

The  format  for  the  geometry  file  follows: 


Zs  Zr  Rhz  Bearl  Bear2  Resolve 


where 

Zs  -  height  of  the  source  in  meters. 

Zr  -  height  of  the  receiver  in  meters. 

Rhz  -  the  horizontal  distance  between  the  source  and  receiver  in 

meters. 

Bearl  -  the  initial  azimuth  for  scanning  in  degrees. 

Bear2  -  the  final  azimuth  for  scanning  in  degrees. 

Resolve  -  the  angular  resolution  of  the  scan  in  degrees.  The  finer  the 
angular  resolution,  the  slower  the  model  will  run.  This  can 
be  overcome  by  using  a  smaller  angle  to  scan. 
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3.2.4  Frequency 


The  frequency  file  contains  the  information  about  the  frequency  to  use  for  the 
model  calculation. 

The  format  for  the  frequency  file  follows: 

Freq  dBsrc  dBback 
where 

Freq  -  frequency  in  Hertz. 

dBsrc  -  level  of  the  source  at  1  meter  in  decibel. 

dBback  -  level  of  the  background  noise  at  receiver  in  decibel. 

If  zero  is  used  for  dBsrc  and  dBback,  the  output  will  be  the  relative  sound 
pressure  level  instead  of  the  sound  pressure  level. 


3.2.5  Ground 

The  ground  file  contains  the  information  to  calculate  the  complex  ground 
impedance.  The  enclosed  table  can  be  used  to  estimate  these  parameters. 


The  format  for  the  ground  file  follows: 


Nground 

Sigmal  Oml  Pnl  Sfl 

D 

Sigma2  Om2  Pn2  Sf2 

where 


Nground 

Sigmal 

Om 

Pnl 


number  of  layers  (1  or  2). 

flow  resistivity  in  cgs  units  of  the  top  layer. 

porosity  of  the  top  layer. 

grain  shape  factor  of  the  top  layer  (normally  use  0.75). 
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Sfl  -  pore  shape  factor  of  the  top  layer  (normally  use  0.875). 

D  -  depth  of  the  top  layer  in  meters. 

Sigma2  -  flow  resistivity  in  cgs  units  of  the  bottom  layer. 

Om2  -  porosity  of  the  bottom  layer. 

Pn2  -  grain  shape  factor  of  the  bottom  layer  (normally  use  0.75). 

Sf2  -  pore  shape  factor  of  the  bottom  layer  (normally  use  0.875). 

Not  all  of  the  entries  are  required.  If  only  one  ground  layer  is  used,  set 
Nground  =  1,  then  D,  Sigma2,  Om2,  Pn2,  and  Sf2  may  be  omitted  from  the 
file.  If  Nground  =  1,  then  values  for  D,  Sigma2,  Om2,  Pn2,  and  Sf2  may  be 
present  in  the  file;  however,  the  program  will  not  read  any  of  these  values. 

3.3  Explanation  of  User  Interface 


A  better  heading  for  this  section  may  be  the  philosophy  of  the  user  interface. 
As  mentioned  at  the  beginning  of  section  3,  the  user  interface  is  designed  to 
allow  easier  model  running  with  numerous  scenarios  and  a  minimum  number 
of  modifications  to  the  input  files.  This  philosophy  will  be  applied  to  the  series 
of  models  being  adapted  by  ARL/BED’s  Acoustics  Team.  The  current  user 
interface  is  being  implemented  on  all  the  flat-earth  non-turbulent  acoustic 
propagation  models  currently  in  the  ARL/BED  repertoire.  As  the  newer 
models  come  out  with  more  complex  interactions,  a  new  interface  will  be 
developed  to  incorporate  the  new  inputs  required  to  execute  the  new  models. 
The  overall  idea  behind  the  user  interfaces  will  be  maintained.  The  normal 
operation  of  the  models  will  require  only  the  basic  input  requirements.  For 
more  advanced  users,  there  are  debug  flags  available  to  provide  for  additional 
input  parameters,  changing  of  default  parameters,  and  output  of  additional 
information  for  the  model.  This  philosophy  provides  a  user  interface  that  is 
easy  to  use,  yet  allows  for  detailed  scientific  studies  to  be  conducted. 

3.4  How  To  Run  SCAFFIP 

SCAFFIP  is  relatively  easy  to  run.  There  is  a  C-shell  script  included  with  the 
source  code  to  compile  the  executable  on  a  Unix  workstation.  To  run  the 
model,  type  sffp  at  the  command  line.  The  program  will  ask  for  the  name  of 
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the  input  files  mentioned  in  section  3.2.  The  printout  labeled  output  shows 
what  will  be  seen  when  the  model  is  run. 

3.5  Format  of  the  Output 

The  output  from  the  model  is  saved  to  a  file  named  scand.  The  format  of  the 
output  file  follows 

Azl 

N1 

R(l)  dB(l) 

R(2)  dB(2) 

R(N1)  dB(Nl) 

Az2 

N2 

R(l)  dB(l) 

R(2)  dB(2) 

R(N2)  dB(N2) 

Az3 

N3 

R(l)  dB(l) 

R(2)  dB(2) 

R(N3)  dB(N3) 


where 

Az  -  azimuth  of  calculation  in  radians. 

N  -  number  of  range  points. 

R(i)  -  range  in  meters. 

dB(i)  -  sound  level  or  relative  sound  level  at  R  in  decibels. 
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Acronyms  and  Abbreviations 


ARL 

Army  Research  Laboratory 

BED 

Battlefield  Environment  Directorate 

CPIO 

Copy  File  Archives  In  and  Out 

FFP 

Fast  Field  Program 

FFT 

Fast  Fourier  Transform 

SCAFFIP 

Scanning  Fast  Field  Program 

TAR 


Tape  Archive 


Appendix  A 

Source  Code  For  SCAFFIP 
General  Version  1.0 


ADMIT.C 


/* 

*  FINDS  ADMITTANCE  OF  ONE  BRANCH  OF  THE  TRANSMISSION  LINE 

*/ 

^include  "complex.h" 

#include  "ffp.h" 

void  admit(admittance,  gp,  last,  layer,  voltage) 

complex  *admittance,  *voltage; 

int  last; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

complex  phase,  ref_l,  ref_2,  unit,  z,  zl,  z2; 
int  i21,  index,  step,  toggle  =  0; 

/* 

*  direction  to  source  from  last  (+1  =  =  up,  -1  =  =  down) 

*/ 

step  =  gp- >  index_source  >  =  last  ?  1  :  -1; 

/* 

*  from  interface  to  section  beyond  interface  (+l=:=up,0=  =  down) 

*/ 

i21  =  (step  +  1)  /  2; 

/* 

*  admittance  of  last  significant  line 
*/ 

Re(*admittance)  =  Re(layer[last].yO); 

Im(*admittance)  =  Im(layer[last].yO); 

/* 

*  step  through  network  to  section  in  front  of  source 
*/ 

Re(unit)  =  1; 

Im(unit)  =  0; 

for  (index  =  last  +  step;  index  <  =  gp- > index_source  &&  step  ==  1 
|  |  index  >  gp- >  index_source  &&  step  ==  -1;  index  +=  step) 

{ 

/* 

*  set  toggle  and  voltage  if  beyond  the  receiver 
*/ 

if  (index  ==  gp- >  index_detector  +  i21) 

{ 
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toggle  =  1; 
Re(*voltage)  =  1; 
Im(*voltage)  =  0; 


} 

electrical  phase  factor  (gamma  =  =  j  *  kz) 

Re(z)  =  0; 

Im(z)  =  -layer[index]  .thickness; 
cmul(z,  z,  layer[index].kz); 
cexp(phase,  z); 

reflection  coefficient  at  the  end  of  section 

csub(zl,  layer[index].yO,  *admittance); 
cadd(z2,  layer[index].yO,  ^admittance); 
cdiv(ref_l,  zl,  z2); 

reflection  coefficient  at  the  beginning  of  section 

cmul(z,  phase,  phase); 
cmul(ref_2,  ref_l,  z); 

total  voltage  at  the  beginning  of  section 

if  (toggle) 

{ 

cadd(zl,  unit,  ref_2); 
cmul(zl,  *voltage,  zl); 
cadd(z2,  unit,  ref_l); 
cmul(z2,  phase,  z2); 
cdiv(*voltage,  zl,  z2); 

} 

admittance  at  the  beginning  of  section 

csub(zl,  unit,  ref_2); 
cmul(zl,  layer[index].yO,  zl); 
cadd(z2,  unit,  ref_2); 
cdiv(*admittance,  zl,  z2); 


admittance  at  the  source 


AIR.C 


/* 

*  AIR  ABSORPTION  ROUTINE 

* 

*  (1)  ANSI  STANDARD  SI. 26-198X 

*  (2)  H.  E.  Bass  et  al, 

*  Absorption  of  Sound  by  the  Atmosphere, 

*  in  "Physical  Acoustics",  1984. 

*/ 

^include  "complex. h" 

/* 

*  reference  ambient  temperature 

*  293.16  degrees  Kelvin  =  =  20  degrees  Celsius 
*/ 

#define  T20  293.15 

double  air_absorption(frequency,  kelvin,  pressure,  relative  humidity) 
double  frequency,  kelvin,  pressure,  relative_humidity; 

{ 

double  alpha,  humidity,  nrf,  orf,  saturation,  square,  temp; 
temp  =  T20  /  kelvin; 

/* 

*  saturation  pressure 

*  reference  2,  page  169,  equation  72 
*/ 

saturation  =  pow(10.0,  8.422  -  10.05916  *  temp  +  5.023  *  loglO(temp)) 

/* 

*  percent  mole  fraction  of  water  vapor 

*  reference  1,  page  19,  equation  DIO 
*/ 

humidity  =  relative_humidity  *  saturation  /  pressure; 

/* 

*  oxygen  relaxation  frequency 

*  reference  1,  page  7,  equation  9 
*/ 

orf  =  pressure  *  (24  +  40400  *  humidity 
*  (humidity  +  0.02)  /  (humidity  +  0.391)); 

/* 

*  nitrogen  relaxation  frequency 

*  reference  1,  page  7,  equation  10 


nrf  =  pressure  *  sqrt(temp) 

*  (9  +  280  *  humidity  *  exp(4.17  *  (1  -  pow(temp,  0.3333)))); 

/* 

*  absorption  coefficient  alpha  (nepers/m) 

*  reference  1,  page  7,  equation  8 
*/ 

square  =  frequency  *  frequency; 

alpha  =  0.01275  *  exp(-2239.1  /  kelvin)  /  (orf  +  square  /  orf); 
alpha  +=  0.1068  *  exp(-3352.0  /  kelvin)  /  (nrf  +  square  /  nrf); 
alpha  =  square 

*  (1.84e-ll  /  sqrt(temp)  /  pressure  +  alpha  *  pow(temp,  2.5)); 
return  alpha; 


CLIP.C 

/* 

*  SIMPLIFY  TRANSMISSION  LINE  BY  CLIPPING  TECHNIQUE  OF  Lee,  et  al 

*/ 

#include  "complex.h" 

#include  "ffp.h" 

/* 

*  exp(kzmax)  =  =  le8 
*/ 

#define  kzmax  18.42068074395236547214 
void  parametersO; 

void  clip(gp,  k,  last,  layer,  load,  njayers,  omega) 

double  k,  omega; 

int  *last,  load,  njayers; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

double  decay  =  0; 

int  i21,  index,  step,  toggle; 

*last  =  -1; 

/* 

*  set  toggle  if  detector  is  not  between  load  and  source 

*/ 

toggle  =  gp->index_detector  >  gp- >  index_source 

?  gp-  >  index_source  >  load  :  gp-  >  index_source  <  =  load; 

/* 

*  direction  to  load  from  source  (  +  1  =  =  up,  -1  =  =  down) 

*/ 

step  =  gp- >  index_source  <  =  load  ?  1  ;  -1; 

/* 

*  conversion  from  interface  to  opposite  layer  (+1  =  =  up,  0  =  =  down) 

*/ 

i21  =  (step  +  1)  /  2; 

/* 

*  step  from  source  to  layer  nearest  load 
*/ 

for  (index  =  gp- >  index_source  +  i21;  index  <  load  &&  step  ==  1 
|  |  index  >  load  &&  step  ==  -1;  index  +=  step) 

{ 

/* 

*  find  admittance  and  wave  number  for  section 
*/ 


parameters(gp,  k,  index,  layer,  njayers,  omega); 

/* 

*  sum  the  exponential  decay  coefficients 
*/ 

decay  -=  layer[index]  .thickness  *  Im(layer[index].kz); 

/* 

*  check  if  beyond  the  detector 
*/ 

if  (toggle) 

{ 

/* 

*  replace  with  infinite  section  if  decay  too  large 
*/ 

if  (decay  >  =  kzmax  /  2) 

{ 

*last  =  index; 
return; 

} 

} 

else 

{ 

/* 

*  v  =  =  0  if  decay  too  large 
*/ 

if  (decay  >  =  kzmax) 
return; 

/* 

*  set  toggle  if  detector  is  reached 

*/ 

if  (index  +  step  ==  gp- >  index_detector  +  i21) 
toggle  =  1; 

} 

} 

/* 

*  find  parameters  for  the  load  layer 
*/ 

parameters(gp,  k,  load,  layer,  njayers,  omega); 

/* 

*  last  significant  layer  in  the  network 
*/ 

*last  =  load; 

} 
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CUTOFF.C 


#include  <stdio.h> 

#include  <stddef.h> 

^include  "complex.h" 

#include  "ffp.h" 

#define  PANEL  VALUES  dflags[5] 
extern  int  ermo; 

void  cutoff(delta_k,  delta_r,  frequency,  gp,  k_max,  k_width,  layer,  n_panels, 
omega,  points,  range,  dflags) 

double  *delta_k,  *delta_r,  frequency,  *k_max,  *k_width,  omega,  range; 
int  *n_panels,  points,  *  dflags; 
struct  gp  *gp; 
struct  layer  *layer; 

{ 

double  c  rnin,  difference  =  0,  gamma; 
int  index,  max,  min,  nyquist; 

FILE  *fp_pv; 


/* 

*  highest  and  lowest  of  layers  between  detector  and  source 
*/ 

max  =  gp-  >  indexdetector  >  gp-  >  index_source 
?  gp-  >  index  detector  :  gp- >  index_source; 
min  =  gp-  >  index  detector  -I-  gp-  >  index_source  -  max  +  1; 

/* 

*  minimum  speed  from  source  to  detector 

*/ 

crnin  =  layer[min].c; 

for  (index  =  min;  index  <  =  max;  index +  +) 

{ 

if  (c_min  >  layer[index].c) 
c_min  =  layer[index].c; 

/* 

*  total  altitude  difference  between  detector  and  source 
*/ 

difference  +=  layerfindex]  .thickness; 

} 

/* 

*  find  upper  cutoff  wave  number  kmax  based  on  an  empirical  relationship 
*/ 


gamma  =  0.0075  *  frequency  +  (2.5  -  6.25e-4  *  frequency)  /  difference; 

*k_max  =  sqrt(gamma  *  gamma  +  (omega  /  c_min)  *  (omega  /  c_min)); 

/* 

*  determine  the  number  of  panels 
*/ 

nyquist  =  *k_max  *  range  /  points  /  PI  +  1; 
if  (*n_panels  =  =  0) 

*n_panels  =  nyquist; 
if  (*n_panels  <  0) 

*n_panels  *=  -nyquist; 

/* 

*  index  of  the  point  nearest  range  based  on  Nyquist  criteria 
*/ 

index  =  *k_max  *  range  /  *n_panels  /  TWO_PI  +  1.5; 

/* 

*  adjust  range  separations  so  that  the  index  point  is  at  range 
*/ 

*delta_r  =  range  /  (index  -  1); 

/* 

*  adjust  upper  cutoff  wave  number  according  to  adjusted  delta_r 
*/ 

*k_max  =  TWO  PI  *  *n_panels  /  *delta_r; 

/* 

*  wave  number  band  width  per  panel 
*/ 

*k_width  =  *k_max  /  *n_panels; 

*delta_k  =  *k_width  /  points; 

/* 

*  Write  Out  Panel  Values 

*/ 

if  (PANELVALUES) 

{ 

if  ((fpjpv  =  fopen("npan.e","a"))  ==  NULL) 

{ 

fprintf(stderr, "error  number  %d\n",ermo); 

perror("npan.e"); 

exit(l); 

} 

fprintf(fp_pv,  "Kmax  =  %10.41f\tNyq  =  %5d\tNpan  =  %5d\n” 

,  *k_max, nyquist ,  *n_panels) ; 

fprintf(fp_pv,"Kwidth  =  %  10.41f,tDelR  =  %10.41f\tDelK  =  %10.41f\n" 
,  *k_width,  *delta_r,  *delta_k) ; 
fclose(fp_pv); 

} 

} 
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FFP.C 

/* 

*  FAST  FIELD  PROGRAM 

*/ 

#include  <stddef.h> 

#include  <stdio.h> 

^include  < stdlib.h  > 

#include  "complex.h" 
include  "ffp.h" 

^define  N  2048 
^define  SQ(X)  (X)  *  (X) 

#define  WAVESPECT  dflags[7] 

#define  SMOOTH  dflags[8] 

void  cutoff(); 
void  fourier(); 
void  profile(); 
void  setup(); 
void  voItage(); 
void  wavenumber(); 
void  zeffectiveO; 
extern  int  ermo; 

double  propmod(a,  fd,  fp,  frequency,  src_level,  back_level,  geometry,  ground, 
n_ground_layers,  n_interfaces,  dflags) 

FILE  *fp; 

double  frequency,  src_level,  back_level; 

int  fd,  n_ground_layers,  n_interfaces,  *dflags; 

struct  geometry  *geometiy; 

struct  ground  *ground; 

struct  interface  *a; 

{ 

complex  amplitude,  field[N],  impedance,  pressure[N],  z,  zO,  zl; 
double  delta_k,  delta_r,  detector,  gain[N],  k,  k_max,  k_min,  k_width, 
omega,  r,  range,  source,  temp,  extra; 
float  buffer[2]; 

int  i,  j,  limit,  n_layers,  n_panels,  points; 
struct  gp  gp; 
struct  layer  *layer; 

FILE  *fp_wspc; 


Preliminary  Setup 


*/ 

detector  =  geometry- >  Zr; 
omega  =  TWO_PI  *  frequency; 
range  =  geometry- >  range; 
source  =  geometry- >Zs; 

/* 

*  Setup  The  Points  and  Panels  For  Transform 
*/ 

setup(&n_panels  ,&points ,  &extra,  dflags) ; 

/* 

*  Impedance  and  Wavenumber  of  Half-Spaces  #1  and  #N 
*/ 

zeffective(a[l].c,  frequency,  ground,  &impedance,  &gp.k_l, 
n_ground_layers,  dflags); 

Re(gp .  impedance_  1 )  =  Re(impedance); 

Im(gp.impedancel)  =  -Im(impedance); 

Im(gp.k_l)  =  -Im(gp.k_l); 

Re(gp.impedancen)  =  Re(gp.k_n)  =  1; 

Im(gp.impedancen)  =  Im(gp.k_n)  =  0; 

/* 

*  Initialize  The  Profile 
*/ 

layer  =  (struct  layer  *)  calloc((size_t)  n_interfaces  +  3, 
sizeof(struct  layer)); 

proftle(a,  detector,  frequency,  &gp,  layer,  n  interfaces,  &n  layers, 
source); 

/* 

*  Squares  of  Intrinsic  Wave  Numbers  Within  Each  Layer 

*/ 

wavenumber(extra,  &gp,  layer,  n_layers,  omega); 

/* 

*  Upper  Cut-Off  Wave  Number 
*/ 

cutoff(&delta_k,  &delta_r,  frequency,  &gp,  &k_max,  &k_width,  layer, 
&n_panels,  omega,  points,  range,  dflags); 

/* 

*  Zero  The  Pressure  Array 
*/ 

for  (i  =  1;  i  <  =  points;  i++) 

Re(pressure[i])  =  Im(pressure[i])  =  0; 

/* 

*  zO  =  (1  -  i)  /  (2  *  sqrt(PI)) 


*/ 

Re(zO)  =  ONE2SQRTPI;  Im(zO)  =  -ONE_2SQRTPI; 

/* 

*  Overlap  And  Add  Field  From  Each  Panel 
*1 

for  (i  =  1;  i  <  =  n_panels;  i++) 

{ 

/* 

*  Starting  Wave  Number  of  Each  Panel 
*/ 

k  =  k_min  =  (i  -  1)  *  k_width; 

/* 

*  Find  The  Pressure  Field  in  The  Transform  Domain 

*/ 

for  (j  =  1;  j  <  =  points;  j  +  +) 

{ 

/* 

*  Pressure  Amplitude 

*/ 

voltagef&amplitude,  &gp,  k,  layer,  njayers,  omega); 
cmul(fteld[j],  amplitude,  zO); 
temp  =  sqrt(k  +  le-20); 

Re(fteld[j])  /-  temp; 

Im(fteld[j])  /=  -temp; 

/* 

*  Write  Out  Horizontal  Wave  Number  Spectrum 

*/ 

if  (WAVESPECT) 

{ 

if  ((fpwspc  =  fopen("wavnum.e","a"))  ==  NULL) 

{ 

fprintf(stderr, "error  number  %d\n", 
ermo); 

perror( "  wavnum .  e " ) ; 
exit(l); 

} 

fprintf(fp_wspc, "  %  10.41f\t%  10.41f\t", 
k,cmod(  amplitude)); 
fprintf(fp_wspc, "  %  10.41f\t%  10.41f\n" , 

Re(amplitude)  ,Im(amplitude)) ; 
fclose(fp_wspc); 


} 

k  +  =  delta_k; 
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} 

/* 

*  Fourier  Transform  The  Wave  Number  Amplitudes 
*/ 

fourier(field,  points,  1); 

/* 

*  Retrieve  The  Incremental  Pressure  Amplitudes 
*/ 

r  =  delta_r; 

for  (j  =  2;  j  <  =  points;  j  +  +) 

{ 

Im(field[j])  =  -Im(field[j]); 

Re(zl)  =  extra  *  sqrt(SQ(r)  +  SQ(source  -  detector)); 
Im(zl)  =  -k_min  *  r; 
cexp(z,  zl); 

cmul(amplitude,  field[j] ,  z); 
temp  =  delta_k  /  sqrt(r); 

Re(amplitude)  *=  temp; 

Im(amplitude)  *=  temp; 

/* 

*  Accumulate  The  Total  Pressure  in  Pressure[j] 

*/ 

Re(pressure[j])  +=  Re(amplitude); 

Im(pressure[j])  +=  Im(amplitude); 
if  (i  =  =  njanels  &&  r  >  range) 
break; 

r  +  =  delta_r; 

} 

} 

limit  =  j  -  2; 
if  (fp  !  =  NULL) 

{ 

fprintf(fp,  "%d\n",  limit); 

} 

else 

{ 

bufferfO]  =  limit; 

write(fd,  (char  *)  buffer,  sizeof(float)); 

} 

limit+  +  ; 
r  =  delta_r; 

for  (i  =  2;  i  <  =  limit;  i++) 

{ 
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*  Find  Gain  Relative  To  0  dB  at  1  meter  (Transmission  Loss) 
*/ 

gain  =  20  *  loglO(cmod(pressure[i])); 

} 

if  (!  SMOOTH) 

{ 

/* 

*  Perform  Smoothing  Algorithm  On  Points 
*/ 

smooth(gain,  limit,  dflags); 

} 

r  =  delta_r; 

for  (i  =  2;  i  <  =  limit;  i++) 

{ 

/* 

*  Write  Out  The  Levels  Versus  Range 

*/ 

gain[i]  =  src_level  +  gain[i]; 
if  ((gain[i]  <  backjevel)  &&  (src_level  >  0)) 
gain[i]  =  backjevel; 
if  (fp  !  =  NULL) 

{ 

fprintf(fp,  "%131f  %131f\n”,  r,  gain[i]); 

} 

else 

{ 

buffer[0]  =  r; 
buffer[l]  =  gain[i] ; 

write(fd,  (char  *)  buffer,  sizeof(double)); 

} 

r  +  =  delta jr; 

} 

return  gain[limit] ; 


} 


FOURIER.C 


/* 

*  Replaces  data  by  its  discrete  Fourier  transform  if  sign  =  =  1  or 

*  by  n  times  its  inverse  discrete  Fourier  transform  if  sign  =  =  -1. 

*  data  is  a  complex  array  of  length  n. 

*  n  must  be  an  integral  power  of  2  (this  is  not  checked  for!). 

* 

*  W.  H.  Press  et  al, 

*  Numerical  Recipes,  Cambridge  University,  1986,  page  394. 

*/ 

#include  "complex.h" 

void  fourier(data,  n,  sign) 
complex  *data; 
int  n,  sign; 

{ 

complex  w,  wp,  ztemp; 

double  temp,  theta; 

int  index,  j  =  1,  m,  max,  step; 


/* 

*  bit-reversal  section 
*/ 

for  (index  =  1;  index  <  =  n;  index++) 

{ 

if  (j  >  index) 

{ 

Re(ztemp)  =  Re(data[j]); 
Im(ztemp)  =  Im(data[j]); 
Re(data[j])  =  Re(data[index]); 
Im(data[j])  =  Im(data[index]); 
Re(data[index])  =  Re(ztemp); 
Im(data[index])  =  Im(ztemp); 

} 

m  =  n  /  2; 

while  (m  >  =  2  &&  m  <  j) 

{ 

j  -=  m; 
m  /=  2; 

} 

j  +  =  m; 
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/* 

*  beginning  of  Danielson-Lanczos  section 
*/ 

max  =  1; 
while  (max  <  n) 

{  /*  while  loop  executed  log2(n)  times  */ 

/* 

*  initialize  for  trigonometric  recurrence 
*/ 

step  =  2  *  max; 

theta  =  PI  /  (max  *  sign); 

temp  =  sin(theta  /  2); 

Re(wp)  =  -2  *  temp  *  temp; 

Im(wp)  =  sin(theta); 

Re(w)  =  1; 

Im(w)  =  0; 

/* 

*  two  nested  inner  loops 

*/ 

for  (m  =  1;  m  <  =  max;  m++) 

{ 

for  (index  =  m;  index  <  =  n;  index  +  =  step) 

{ 

/* 

*  Danielson-Lanczos  formula 
*/ 

j  =  index  +  max; 
cmul(ztemp,  data[j],  w); 
csub(data[j],  data[index],  ztemp); 
cadd(data[index] ,  data[index],  ztemp); 

} 

/* 

*  trigonometric  recurrence 
*/ 

cmul (ztemp,  w,  wp); 
cadd(w,  w,  ztemp); 

} 

max  =  step; 

} 

} 
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HEADER. C 

^include  "ffp.h" 

#include  "complex.h" 

/* 

*  print  a  line  of  80  dashes 
*/ 

#defme  dash\ 


printff - 

printff" - \n") 


#define  RADTODEG(x)  (180.  *  x)  /  PI 

void  header(a,  frequency,  geometry,  ground,  n_ground_layers,  n_interfaces) 

double  frequency; 

int  n_ground  layers,  n_interfaces; 

struct  geometry  ^geometry; 

struct  ground  Aground; 

struct  interface  *a; 

{ 

int  index; 

printff "\n\n\n  SINGLE  FREQUENCIES\n”); 

dash; 

printff  "RECEIVER  HEIGHT  (Zr)  =  %7.  Ilf  METERS  SOURCE  HEIGHT  (Zs)  =  %8.11f 
METERSXn",  geometry- >  Zr,  geometry- >Zs); 

printff'RANGE  =  %7.  Ilf  METERS  BEARING  OF  SOURCE  =%lf 

DEGREES\n" , 

geometry-  >  range,  RAD_TO_DEG(geometry-  >  theta_  1 )) ; 
printff'FREQUENCY  =  %7d  Hz\n”,  (int)  frequency); 

dash; 

printf("  KELVIN  RELATIVE  ATMOSPHERIC  WIND 

WIND\n"); 

printff  "INTERFACE  HEIGHT  TEMPERATURE  HUMIDITY  PRESSURE  SPEED 
DIRECTION^"); 

for  (index  =  n_interfaces;  index  >  0;  index-) 

printff"  %6d%  13.21f%  13.21f%  12.  llf%  13.31f%  10.31f%  1 1 .31f\n\ 
index,  a[index].z,  a[index] .kelvin, 
a[index]  ,relative_humidity ,  a[index]  .pressure, 
a[index] .  wind_speed,  a[index] .  wind_theta); 

dash; 

if  (n_ground_layers  =  =  1) 

{ 
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printf("  GROUND  LAYER  SIGMA  POROSITY\n"); 
printf("  l%15.21f%13.21f\n", 

ground-  >  sigma_  1 ,  ground-  >  omega_  1 ) ; 


{ 

printf(" GROUND  LAYER  SIGMA  POROSITY  DEPTHVn"); 

printf("  1  %  15.21f%  13.21f%  10.31f\n", 

ground- >sigma_l,  ground- > omega_l ,  ground- > depth); 
printf("  2%15.21f%13.21Tn", 

ground-  >  sigma_2 ,  ground-  >  omega_2) ; 

} 

dash; 

} 
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INPUT. C 


^include  <stddef.h> 

^include  <stdio.h> 

^include  < stdlib.h > 

^include  "ffp.h" 

^include  "complex.h" 

#defme  MET  FORM  dflags[l] 

#defme  openfile(NAME)\ 

printf("Input  name  of  NAME  file:\t");\ 
scanf("%s",  filename);\ 
printf("%s\n",  filename);\ 
if  ((fp  =  fopen(filename,  "r"))  =  =  NULL)  {\ 
fprintf(stderr,  "error  number  %d\n",  ermo);\ 
perror(filename);\ 
exit(l);\ 

} 

^define  DEG  TO  RAD(x)  (PI  *  x)  /  180. 
extern  int  ermo; 

void  input(a,  dflags,  frequency,  srcjevel,  backjevel,  geometry,  ground, 
n  groundjay ers ,  ninterfaces ,  n_freq) 
int  *n_ground_layers,  *n_interfaces,  *n_freq,  *dflags; 
double  *  frequency,  *src_level,  *back_level; 
struct  geometry  *geometry; 
struct  ground  Aground; 
struct  interface  **a; 

{ 

FILE  *fp; 
char  filename[81]; 
int  code, index; 

double  wind_theta,azmuth  1  ,azmuth2  .dazmuth; 

/* 

*  Read  Debug  File 

*/ 

openfile(debug); 
fscanf(fp,  "%d",  &code); 
fclose(fp); 
parser(code ,  dflags) ; 

/* 

*  Read  Weather  File 
*/ 


openfile(weather); 
fscanf(fp,  "%d",  n_interfaces); 

*a  =  (struct  interface  *)  ca!loc((size_t)  *n_interfaces  +  1 , 
sizeof(struct  interface)); 
if  (METFORM) 

{ 

/* 

*  Wind  Direction  in  Radians 

*/ 

for  (index  =  1;  index  <  =  *n_interfaces;  index ++) 

{ 

fscanf(fp,  "%lf  %lf  %lf  %lf  %lf  %lf", 

&(*a)  [index]  .z,  &(*a)[index]  .kelvin, 

&(*a)[index]  .relative_humidity , 

&(*a)[index]  .pressure,  &(*a)[index]  ,wind_speed, 
&(*a)[index]  ,wind_theta); 

} 

} 

else 

{ 

/* 

*  Wind  Direction  in  Degrees 
*/ 

for  (index  =  1;  index  <=  *n_interfaces;  index ++) 

{ 

fscanf(fp,  "%lf  %lf  %lf  %lf  %lf  %lf”, 

&(*a)[index]  .z,  &(*a)[index] .kelvin, 

&(*a)[index]  ,relative_humidity , 

&(*a)[index]  .pressure, 

&(*a)[index]  ,wind_speed,  &wind_theta); 
(*a)[index].wind_theta  =  DEG_TO_RAD(wind_theta); 

} 

} 

fclose(fp); 

/* 

*  Read  Geometry  File 
*/ 

openfile(geometry) ; 

fscanf(fp,  "%lf  %lf  %lf  %lf  %lf  %lf". 

&(geometry-  >  Zs),  &(geometry->Zr),  &(geometry->  range), 
&azmuthl,  &azmuth2,  &dazmuth); 
geometry- >  theta_l  =  DEG_TO_RAD(azmuthl); 
geometry- >  theta_2  =  DEG_TO_RAD(azmuth2); 
geometry- >  delta_theta  =  DEG_TO_RAD(dazmuth); 


fclose(fp); 


/* 

*  Read  Frequency  File 
*/ 

openfile(frequency ) ; 

/* 

*  Srcjevel  and  Backjevel  is  local  to  input. c 
*/ 

fscanf(fp,  "%lf  %lf  %lf",  frequency,  src_level,  backjevel); 
fclose(fp); 

/* 

*  Read  Ground  File 
*/ 

openfile(ground); 

fscanf(fp,  "  %d  %lf  %lf  %lf  %lf", 

n  ground  layers,  &(ground->sigma_l),  &(ground->omega_l), 
&(ground-  >  pn_l),&(ground-  >  sf_l)); 
if  (*n_ground  layers  >  1) 

fscanfffp,  "%lf  %lf  %lf  %lf  %lf",  &(ground->  depth), 
&(ground-  >  sigma_2),  &(ground-  > omega_2), 

&(ground-  >  pn_2),&(ground-  >  sf_2)); 
fclose(fp); 


INSERT. C 


/* 

*  insert  a  new  layer  with  top  interface  height  =  =  height 
*/ 

^include  "ffp.h" 

void  insert(height,  index,  layer,  njayers,  z) 
double  height,  *z; 
int  *  index,  *n_layers; 
struct  layer  *layer; 

{ 

/* 

*  find  the  index  of  height  in  z  if  it  exists 
*/ 

for  (*index  =  1;  *index  <  *n_layers;  +  +  *index) 
if  (z[*indexj  =  =  height) 

/* 

*  return  the  index  of  z  that  yields  height 


*  shift  array  z  up  to  open  a  slot  for  height 
*/ 

for  (*index  =  *n_layers;  *index  >  1;  ~*index) 

{ 

layer[*index+l].c  =  layer[*index].c; 
layer[*index+l].mu  =  layer[*index].mu; 
layer[*index+l].rho  =  layer[*index].rho; 
z[*index+l]  =  z[*index]; 

/* 

*  check  if  height  belongs  above  the  next  z 
*/ 

if  (height  >  z[*index-l]) 
break; 

} 

/* 

*  extend  the  array  length  and  insert  height 
*/ 

+  +*n_layers; 
z[*index]  =  height; 


} 


MAIN.C 


/* 

*  M ASTER. C  -  Main  Driver  For  General  Purpose  User  Interface 

* 

*  John  M.  Noble 

*  U.S.  Army  Research  Laboratory 

*  Battlefield  Environment  Directorate 

*  ATTN:  AMSRL-BE-M 

*  White  Sands  Missile  Range,  New  Mexico 

* 

*  History: 

* 

*  Dec  1991  -  Converted  to  C  by  Frank  Owens  -  Ball  State  Univ. 

*  Feb  1993  -  User  Interface  Generalized 

*/ 

^include  <fcntl.h> 

^include  <stdio.h> 

^include  <  stdlib.h > 

^include  "complex.h" 

^include  "ffp.h" 

^define  FLAGS  O  WRONLY  |  O  CREAT  |  0_TRUNC 
#define  MODE  0644  /*  this  is  the  file  mode  for  a  created  scan.bin  */ 

^define  RAD  TO  DEG(x)  (180.  *  x)  /  PI 

extern  int  ermo; 

double  propmod(); 
void  header(); 
void  input(); 
void  parser(); 
void  sound(); 

main(argc,  argv) 
int  argc; 
char  **argv; 

{ 

FILE  *fp  =  NULL; 

char  *id  =  ”@(#)  December  1991  version  -  Frank  W.  Owens”; 
double  attenuation,  frequency,  src  level,  back_level,  theta; 
float  buffer; 

int  fd,  index,  n_ground_layers,  n_interfaces,  dflags[15]; 
struct  geometry  geometry; 


struct  ground  ground; 
struct  interface  *a; 
if  (argc  =  =  1) 

{ 

if  ((fp  =  fopen("scand",  "w"))  —  NULL) 

{ 

fprintf(stderr,  "error  number  %d\n",  ermo); 

perror("scan"); 

exit(l); 

} 

} 

else  if  (argc  =  =  2) 

{  . . 

if  ((fd  =  open("scand.bin",  FLAGS  ,  MODE))  <  0) 

{ 

fprintf(stderr,  "error  number  %d\n",  ermo); 

perror(  "scan.bin"); 

exit(l); 

} 

} 

else 

{ 

fprintf(stderr,  "usage:  %s  [b]\n”,  *argv); 
exit(2); 

} 

input(&a,  dflags,  &frequency,  &src_level,  &back_level,  &geometry, 
&ground,  &n_ground_layers,  &n_interfaces); 
header(a,  frequency,  &geometry,  &ground,  n_ground_layers,  n_interfaces), 

/* 

*  Begin  Scanning  Through  The  Azimuths 
*/ 

for  (theta  =  geometry. theta_l;  theta  <  =  geometry ,theta_2; 
theta  +=  geometry.  delta_theta) 

{ 

if  (argc  =  =  1) 

{ 

fprintf(fp,  "%lf\n",  theta); 

} 

else 

{ 

buffer  =  theta; 

write(fd,  (char  *)  &buffer,  sizeof(float)); 

} 

sound(a,  n_interfaces,  theta,  dflags); 


attenuation  =  propmod(a,  fd,  fp,  frequency,  srclevel, 
back_level,  &geometry,  &ground,  n_ground_layers, 
n_interfaces,  dflags); 

printf(" Position  Angle  =  %lf\n",  RAD_TO_DEG(theta)); 
printf("%8.21f  %9.41f\n",  geometry. range,  attenuation); 


} 

if  (argc  =  =  1) 
fclose(fp); 

else 

close(fd); 

printf(" - ”); 

printfC - \n"); 

return  0; 


PARAMETE.C 

- 

*  NETWORK  TRANSMISSION  LINE  PARAMETERS  FOR  SECTION  l  index 

*/ 

^include  "complex.h"  . 

^include  "ffp.h" 

void  parameters(gp,  k,  l_index,  layer,  n_layers,  omega) 

double  k,  omega; 

int  l_index,  n_layers; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

complex  z; 

/* 

*  electrical  wave  number  kz  =  =  -  j  *  gamma 
*/ 

Re(z)  =  Re(layer[l_index].ki2)  -  k  *  k; 

Im(z)  =  Im(layer[l_index].ki2); 
csqrt(layer[l_index].kz,  z); 

/* 

*  Re(gamma)  is  positive  for  the  correct  branch 
*/ 

if  (Im(layer[l_index].kz)  >  0) 

{ 

Re(layer[l  index], kz)  =  -Re(layer[l_index].kz); 

Im(layer[l_index].kz)  =  -Im(layer[l_index].kz); 

} 

/* 

*  electrical  characteristic  admittance 
*/ 

Re(layer[l_index]  ,y0)  = 

Re(layer[l_index].kz)  /  (omega  *  layer[l_index].rho); 

Im(layer[l_index].yO)  = 

Im(layer[l_index].kz)  /  (omega  *  layer[!_index].rho); 

/+ 

*  admittance  for  half-spaces 
*/ 

if  (l_index  =  =  1) 

{ 

cmul(z,  gp-  >  impedance_l ,  gp->k_l); 
cdiv(layer[l_index].yO,  layer[l_index].yO,  z); 
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} 

if  (l_index  =  =  n_layers) 

{ 

cmul(z,  gp-  >  impedance_n,  gp->k_n); 
cdiv(layer[l_index].yO,  layer[l_index].yO,  z); 


} 


PARSER.  C 


^include  <stdio.h> 


void  parser(code,dflags) 
int  code,  *dflags; 

{ 

/* 

*  DFLAGS(O)  -  2A0  =  1 

*/ 

dflagsfO]  =  (1  &  code); 

/* 

*  DFLAGS(l)  -  2A1  =  2 

*/ 

dflags[l]  =  (2  &  code)  /  2; 

/* 

*  DFLAGS(2)  -  2A2  =  4 

*/ 

dflags[2]  =  (4  &  code)  /  4; 

/* 

*  DFLAGS(3)  -  2A3  =  8 
*/ 

dflags[3]  =  (8  &  code)  /  8; 

/* 

*  DFLAGS(4)  -  2A4  =  16 
*/ 

dflags[4]  =  (16  &  code)  /  16; 

/* 

*  DFLAGS(5)  -  2A5  =  32 
*/ 

dflags[5]  =  (32  &  code)  /  32; 

/* 

*  DFLAGS(6)  -  2A6  =  64 
*/ 

dflags[6]  =  (64  &  code)  /  64; 

/* 

*  DFLAGS(7)  -  2A7  =  128 

*/ 

dflags[7]  =  (128  &  code)  /  128; 

/* 

*  DFLAGS(8)  -  2A8  =  256 
*/ 

dflags[8]  =  (256  &  code)  /  256; 
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*  DFLAGS(9)  -  2A9  =  512 
*/ 

dflags[9]  =  (512  &  code)  /  512; 

/* 

*  DFLAGS(IO)  -  2A10  =  1024 
*/ 

dflags[10]  =  (1024  &  code)  /  1024; 

/* 

*  DFLAGS(ll)  -  2A11  =  2048 

*/ 

dflagsfll]  =  (2048  &  code)  /  2048; 

/* 

*  DFLAGS(12)  -  2A12  =  4096 

*/ 

dflags[12]  =  (4096  &  code)  /  4096; 

/* 

*  DFLAGS(13)  -  2A13  =  8192 

*/ 

dflags[13]  =  (8192  &  code)  /  8192; 

/* 

*  DFLAGS(14)  -  2A14  =  16384 
*/ 

dflags[14]  =  (16384  &  code)  /  16384; 


} 


PROFILE.C 

/* 

*  INITIALIZE  PROBLEM  ATMOSPHERIC  PROFILE 

*/ 

#include  <stddef.h> 

#include  < stdlib.h > 

^include  "ffp.h" 

double  air_absorption(); 
void  insert(); 


void  profile(a,  detector,  frequency,  gp,  layer,  n_interfaces,  njayers,  source) 

double  detector,  frequency,  source; 

int  n_interfaces,  *n_layers; 

struct  interface  *a; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

double  *z; 
int  index; 

/* 

*  initialize  the  problem  profile  (z  increases) 

*/ 

z  =  (double  *)  calloc((size_t)  n_interfaces  +  3,  sizeof(double)); 
for  (index  =  1;  index  <=  n_interfaces;  index ++) 

{ 

layer[index].c  =  a[index].c; 

layer[index].mu  =  air_absorption(frequency,  a[index].kelvin, 
afindex]  .pressure,  a[index]  ,relative_humidity); 
layer[index].rho  =  1.2; 
z[index]  =  a[index].z; 

} 

*njayers  =  n_interfaces; 

/* 

*  insert  detector  and  source  interfaces  (lowest  first) 

*/ 

if  (source  <  detector) 

{ 

insert(source,  &gp->index_source,  layer,  njayers,  z); 
insert(detector,  &gp- >  index_detector,  layer,  njayers,  z); 

} 

else 
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{ 


insert(detector,  &gp- >  index_detector,  layer,  njayers,  z); 
insert(source,  &gp- >  index_source,  layer,  njayers,  z); 


convert  interface  heights  to  layer  thicknesses 

for  (index  =  2;  index  <  *njayers;  index ++) 
layer[index]  .thickness  =  z[index]  -  z[index-l]; 


SETUP.C 

#include  <stddef.h> 

^include  <stdio.h> 

#include  < stdlib.h > 

#defme  DEFAULTEXTRALOSS  le-4 
#define  EXTRALOSSFLAG  dflags[3] 

#define  POINTS  1024 

#define  N_PANELS  -2 

^define  TRANSFORMFLAG  dflags[4] 

extern  int  ermo; 

void  setup(n__panels, points, extra, dflags) 
double  *extra; 

int  *n_panels,  *points,  *dflags; 

{ 

FILE  *fp; 
char  filename[81]; 

/* 

*  Extra  Loss 
*/ 

if  (EXTRA  LOSS  FLAG) 

{ 

if  ((fp  =  fopen("extra.los",  ”r"))  ==  NULL) 

{ 

fprintf(stderr,  "error  number  %d\n",  ermo); 

perror("extra.los"); 

exit(l); 

} 

fscanf(fp,  ”%lf",  extra); 
fclose(fp); 

} 

else 

*extra  =  DEFAULTEXTRALOSS; 

/* 

*  Number  of  Points  and  Panels 
*/ 

if  (TRANSFORM  FLAG) 

{ 

if  ((fp  =  fopenC’trans.frm",  ”r"))  ==  NULL) 

{ 

fprintf(stderr,  "error  number  %d\n\  ermo); 


perror( "  trans .  frm " ) ; 
exit(l); 

} 

fscanf(fp,  ”%d  %d",  points, n_panels); 
fclose(fp); 


{ 

*points  =  POINTS; 
*n_panels  =  N_PANELS; 

} 

} 
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SMOOTH.  C 

^include  <stdio.h> 

^define  N  1024 

void  smooth(gain, limit, dflags) 
double  *gain; 
int  limit,  *dflags; 

{ 

int  loop,  i,  npt; 
double  tmp[N],  avg; 

/* 

*  Copy  Data  to  Working  Array 
*/ 

for  (loop  =  0;  loop  <  =  limit;  loop  +  +) 
tmp[loop]  =  gain[loop]; 

/* 

*  Smooth  The  First  Five  Points 

*/ 

npt  =  3; 

for  (loop  =  3;  loop  <  7;  loop++) 

{ 

avg  =  0.; 

for  (i  =  -(npt-l)/2;  i  <  =  (npt-l)/2;  i  +  +) 
avg  +=  tmp[loop+i]; 

/* 

*  Replace  Gain  at  R[loop]  By  3,  5,  7  or  9  Point  Average 
*/ 

gain[loop]  =  avg  /  npt; 
npt  +=2; 

} 

/* 

*  1 1  Point  Average  For  Smoothing 
*/ 

for  (loop  =  7;  loop  <  =  limit-5;  loop  +  +) 

{ 

avg  =  0.; 

for  (i  =  -5;  i  <=5;  i++) 
avg  +=  tmp[loop+i]; 

/* 

*  Replace  Gain  at  R[loop]  By  1 1  Point  Average 
*/ 

gain[loop]  =  avg  /  11.; 


} 


/* 

*  Smooth  The  Last  Five  Points 
*/ 

npt  =  9; 

for  (loop  =  limit-4;  loop  <  =  limit;  loop++) 

{ 

avg  =  0.; 

for  (i  =  -(npt-l)/2;  i  <  =  (npt-l)/2;  i  +  +) 
avg  +=  tmp[loop+i]; 

/* 

*  Replace  Gain  at  R[loop]  By  3,  5,  7  or  9  Point  Average 
*/ 

gain[loop]  =  avg  /  npt; 
npt  -=2; 

} 

} 
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SOUND. C 


/* 

*  calculate  the  sound  speed  profile 
*/ 

^include  <stdio.h> 

#include  "complex.h" 

#include  "ffp.h" 

#deftne  TO  273.16  /*  Kelvin  freezing  point  of  H20  */ 

^define  RO  8314.16  /*  universal  gas  constant  * 

^define  SOUND_SPEED_OUT_FLAGdflags[0] 

^define  RECIPROCITYFLAG  dflags[2] 

extern  int  ermo; 

void  sound(a,  n_interfaces,  bearing,  dflags) 
double  bearing; 
int  n_interfaces,  *dflags; 
struct  interface  *a; 

{ 

FILE  *fp_sound; 
int  index; 

double  gamma,  humidity,  lpsat,  R,  temp,  theta; 

for  (index  =  1;  index  <=  n_interfaces;  index ++) 

{ 

/* 

*  Angle  Between  Wind  And  Propagation  Path 
*/ 

if  (RECIPROCITY  FLAG) 

{ 

/* 

*  Source  Centered  Geometry 
*/ 

theta  =  a[index]  ,wind_theta  -  PI  -  bearing; 

} 

else 

{ 

/* 

*  Sensor  Centered  Geometry 
*/ 

theta  =  a[index]  .wind_theta  -  bearing; 

} 

/* 

*  The  Partial  Pressure  of  Saturated  Water  Vapor 
*/ 

temp  =  a[index].kelvin  /  TO; 


/* 

* 

*/ 


/* 

* 

*/ 


lpsat  =  -10.79586  /  temp  -  5.02808  *  loglO(temp) 

-  1.50474e-4  *  pow(10.0,  8.29692  *  (1  -  temp)) 
4-  4.28730e-4  *  pow(10.0,  4.76955  *  (1  -  temp)) 
+  8.575983444; 

The  Fraction  of  Water  Molecules  in  Air 

humidity  =  a[index].relative_humidity 

*  pow(10.0,  (lpsat  -  2))  /  afindex], pressure; 

Ratio  of  Specific  Heats 

gamma  =  (humidity  +  7)  /  (humidity  4-  5); 

R  =  RO  /  (29  -  11  *  humidity); 


/* 

*  Sound  Speed  Due  to  Humidity  And  Temperature 

*/ 

a[index].c  =  sqrt(a[index].kelvin  *  gamma  *  R); 

/* 

*  Adjustment  For  Wind 

*/ 

a[index].c  4-=  a[index].wind_speed  *  cos(theta); 

} 

/* 

*  Print  Sound  Speed  Profile  If  Desired 

*/ 

if  (SOUNDSPEEDOUTFLAG) 

{ 

if  ((fpsound  =  fopen("sound.o","a"))  =  =  NULL) 

{ 

fprintf(stderr, "error  number  %d\n",ermo); 

perror(''sound.o"); 

exit(l); 

} 

for  (index  =  1;  index  <  =  n_interfaces;  index  4- 4-) 
fprintf(fp_sound, "  %  10.41f\t%  10.41f\n’\ 
a[index]  ,z,a[index]  .c); 
fclose(fp_sound); 

} 

} 
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VOLTAGE. C 


/* 

*  FINDS  THE  VOLTAGE  detector_voltage  AT  index_detector 

*  DUE  TO  CURRENT  SOURCE  source_current  AT  index_source 

*/ 

^include  "complex.h" 

#include  "ffp.h" 

void  admit(); 
void  clip(); 

void  voltage(detector_voltage ,  gp,  k,  layer,  n_layers,  omega) 

complex  *detector_voltage; 

double  k,  omega; 

int  njayers; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

complex  admittancel,  admittance_n,  source_admittance, 
sourcecurrent,  sourcevoltage,  v; 
int  last  1 ,  last  n; 


/* 

*  zero  default  voltage 
*/ 

Re(*detector_voltage)  =  Im(*detector_voltage)  =  Re(v)  =  Im(v) 

/* 

*  find  last  significant  line  toward  half-space  #N 

*/ 

clipfgp,  k,  &last_n,  layer,  njayers,  njayers,  omega); 
if  (last_n  !  =  -1) 

{  I*  if  (v  =  =  0)  return  zero  volts  */ 

/* 

*  find  last  significant  line  toward  half-space  #1 
*/ 

clipfgp,  k,  &last_l,  layer,  1,  njayers,  omega); 
if  (last  J  !  =  -1) 

{  /*  if  (v  =  =  0)  return  zero  volts  */ 

/* 

*  admittance  that  the  source  sees  looking  toward  #N 
*/ 

admit(&admittance_n,  gp,  last_n,  layer,  &v); 

/* 


admittance  that  the  source  sees  looking  toward  #1 

admit(&admittance_l,  gp,  last_l,  layer,  &v); 

total  source  admittance  (add  in  parallel) 

cadd(source_admittance,  admittance_l ,  admittance_n); 

source  current  (equation  14) 

Re(source_current)  =  2  *  k  / 

(layer[gp->index_source].rho  *  omega); 
Im(source_current)  =  0; 

source  voltage  (equation  20) 

cdiv(source_voltage,  source_current,  source_admittance); 
detector  voltage  (equation  23) 
cdiv(*detector_voltage,  source_voltage,  v); 


WAVENUMB.C 


*  SQUARES  OF  COMPLEX  INTRINSIC  WAVENUMBERS  FOR  EACH  LINE  SECTION 

*/ 

^include  "complex.h" 

#include  "ffp.h" 

void  wavenumber(extra,  gp,  layer,  n  layers,  omega) 

double  extra,  omega; 

int  njayers; 

struct  gp  *gp; 

struct  layer  *layer; 

{ 

complex  z; 
int  index; 

/* 

*  finite  thickness  slabs 

*/ 

for  (index  =  2;  index  <  njayers;  index++) 

{ 

Re(z)  =  omega  /  layer[index].c; 

Im(z)  =  -  extra  -  layer[index]  .mu; 
cmul(layer[index].ki2,  z,  z); 

} 

/* 

*  half-spaces  #1  and  #N 
*/ 

Re(z)  =  omega  /  layer[l].c  *  Re(gp->k_l); 

Im(z)  =  omega  /  layer[l].c  *  Im(gp->k_l)  -  extra  -  layer[l].mu; 
cmul(layer[l].ki2,  z,  z); 

Re(z)  =  omega  /  layer[njayers]  .c  *  Re(gp->k_n); 

Im(z)  =  omega  /  layer[njayers].c  *  Im(gp->k_n) 

-  extra  -  layer[njayers].mu; 
cmul(layer[njayers].ki2,  z,  z); 

} 
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Z2LAYERS.C 


/* 

*  FIND  THE  IMPEDANCE  FOR  ONE  OR  TWO  GROUND  LAYERS 

* 

*  K.  Attenborough, 

*  Acoustical  Impedance  Models  for  Outdoor  Ground  Surfaces, 

*  J.  Sound  Vib.  99  (1985),  521-544. 

* 

*  See  equation  15  on  page  527  of  the  above  reference. 

*/ 

#include  "complex.h" 

void  z21ayers(depth,  impedance,  kl,  zbottom,  ztop) 
complex  *impedance,  *kl,  zbottom,  ztop; 
double  depth; 

{ 

complex  factor,  z,  zl,  z2,  ztan; 
double  a,  b; 

/* 

*  Zb  =  zbottom 

*  Zt  =  ztop 

*  Zb  -  i  Zt  tan(a  +  ib)  zl 

*  impedance  = - Zt  =  —  Zt  =  factor  Zt 

*  Zt  -  i  Zb  tan(a  +  ib)  z2 

*/ 

a  =  depth  *  Re(*kl); 
b  =  depth  *  Im(*kl); 

/* 

*  tan(a  +  ib)  =  (tan(a)  +  i  tanh(b))  /  (1  -  i  tan(a)tanh(b)) 

*/ 

Re(zl)  =  tan(a); 

Im(zl)  =  tanh(b); 

Re(z2)  =  1; 

Im(z2)  =  -Re(zl)  *  Im(zl); 
cdiv(ztan,  zl,  z2); 

/* 

*  determine  the  factor 
*/ 

cmul(z,  ztop,  ztan); 

Re(zl)  =  Re(zbottom)  4-  Im(z); 

Im(zl)  =  Im(zbottom)  -  Re(z); 
cmul(z,  zbottom,  ztan); 

Re(z2)  =  Re(ztop)  +  Im(z); 

Im(z2)  =  Im(ztop)  -  Re(z); 
cdiv(factor,  zl,  z2); 

/* 
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complex  impedance 
cmul(*impedance,  factor,  ztop); 


Z4PARAME.C 


/* 

*  FOUR-PARAMETER  APPROXIMATE  MODEL  FOR  GROUND  IMPEDANCE 

* 

*  K.  Attenborough, 

*  Acoustical  Impedance  Models  for  Outdoor  Ground  Surfaces, 

*  J.  Sound  Vib.  99  (1985),  521-544. 

*/ 

^include  "complex.h" 

^define  GAMMA  1.4 

^define  RHOO  1.21  /*  k  /  mA3  */ 

void  z4parameter(c,  frequency,  impedance,  wavenumber,  omega,  pn,  sf,  sigma) 
complex  *  impedance,  *wavenumber; 
double  c,  frequency,  omega,  pn,  sf,  sigma; 

{ 

complex  kb,  ztemp; 
double  alpha,  beta,  q2,  x; 

alpha  =  TWO  PI  *  frequency; 
beta  =  alpha  /  c; 

/* 

*  1000  converts  from  cgs  to  mks 
*/ 

x  =  1000  *  sf  *  sf  *  sigma  /  (RHOO  *  alpha); 

/* 

*  q2  defined  after  equation  8  on  page  524 
*/ 

q2  =  pow(omega,  -pn); 

/* 

*  kb  defined  by  equation  1 1  on  page  524 
*/ 

Re(ztemp)  =  1.5764  *  q2; 

Im(ztemp)  =  GAMMA  *  omega  *  x; 
csqrt(kb,  ztemp); 

/* 

*  impedance  defined  by  equation  12  on  page  524 
*/ 

Re(ztemp)  =  (4  *  q2)  /  (3  *  omega); 

Im(ztemp)  =  x; 
cdiv(*impedance,  ztemp,  kb); 


/* 

*  kb  has  no  units 
*/ 

Re(*wavenumber)  =  beta  *  Re(kb); 
Im(*wavenumber)  =  beta  *  Im(kb); 

} 
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ZEFFECTI.C 


/* 

*  CALCULATES  EFFECTIVE  GROUND  IMPEDANCE 
*/ 

#include  <stdio.h> 

#include  <stddef.h> 

#include  "complex. h" 

^include  "ffp.h” 

^define  PRT  IMPED  dflags[6] 

void  z21ayers(); 
void  z4parameter(); 
extern  int  ermo; 

void  zeffective(c,  frequency,  ground,  impedance,  kl,  n_ground_layers,  dflags) 

complex  *impedance,  *kl; 

double  c,  frequency; 

int  n  ground  layers,  *dflags; 

struct  ground  Aground; 

{ 

complex  k2,  zbottom,  ztop; 

FILE  *fp_imp; 


/* 

*  find  impedance  for  top  ground  interface 
*/ 

z4parameter(c,  frequency,  &ztop,  kl,  ground- >omega_l,  ground->pn_l,  ground- > sf_l , 
ground-  >  sigma_l); 

Re(*impedance)  =  Re(ztop); 

Im(*impedance)  =  Im(ztop); 
if  (n_ground_layers  =  =  2) 

{ 

I* 

*  put  in  the  contribution  from  the  bottom  ground  interface 

*1 

z4parameter(c,  frequency,  &zbottom,  &k2,  ground- >omega_2,  ground->pn_  , 

ground-  >  sf_2, 

ground-  >  sigma_2); 

z21ayers(ground-> depth,  impedance,  kl,  zbottom,  ztop); 

} 
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Write  Out  Impedance  Value 


if  (PRTJMPED) 

{ 

if  ((fp_imp  =  fopen("imped.e","a"))  ==  NULL) 

{ 

fprintffstderr,  "error  number  %d\n\ermo); 

perror("imped.e"); 

exit(l); 

} 

fprintf(fp_imp,  "Impedance  =  (%  10.41f,  %  10.41f)\n” , 
Re(*impedance) ,  -Im(*impedance)) ; 
fprintf(fp_imp,  "Wavenumber  =  (%10.41f,%10.41f)\n", 
Re(*kl),-Im(*kl)); 
fclose(fp_imp); 


COMPLEX.H 

#include  <math.h> 

^define  ONE_2SQRTPI  0.28209479177387814347 
#define  PI  3.14159265358979323846 
#define  TWO_PI  6.28318530717958647692 

tfifndef  _COMPLEX_ 

#define  _COMPLEX_ 

typedef  struct 

{ 

double  real,  imaginary; 

}  complex; 

/* 

*  Re(z)  is  the  real  part  of  z 

*  Im(z)  is  the  imaginary  part  of  z 
*/ 

^define  Re(z)  (z).real 
#define  Im(z)  (z). imaginary 

/* 

*  z  is  the  complex  sum  of  zl  and  z2 
*/ 

^define  cadd(z,  zl,  z2)\ 

Re(z)  =  Re(zl)  +  Re(z2);\ 

Im(z)  =  Im(zl)  +  Im(z2) 

/* 

*  z  is  the  complex  quotient  of  zl  and  z2 

*  z  =  zl  /  z2 
*/ 

#define  cdiv(z,  zl,  z2)\ 

double  denominator,  tempi,  temp2;\ 

\ 

denominator  =  Re(z2)  *  Re(z2)  4-  Im(z2)  *  Im(z2);\ 
tempi  =  Im(zl)  *  Re(z2);\ 
temp2  =  Re(zl)  *  Im(z2);\ 

Re(z)  =  ((Re(zl)  +  Im(zl))  *  (Re(z2)  +  Im(z2))\ 

-  tempi  -  temp2)  /  denominator^ 

Im(z)  =  (tempi  -  temp2)  /  denominator^ 


} 

/* 

*  z  is  the  complex  exponential  of  zl 
*/ 

^define  cexp(z,  zl)\ 

{\ 

double  r,  theta;\ 

\ 

r  =  exp(Re(zl));\ 

theta  =  fmod(Im(zl),  TWO_PI);\ 

Re(z)  =  r  *  cos(theta);\ 

Im(z)  =  r  *  sin(theta);\ 

} 

/* 

*  cmod(z)  returns  | z| 

*/ 

#define  cmod(z)  sqrt(Re(z)  *  Re(z)  +  Im(z)  *  Im(z)) 

/* 

*  z  is  the  complex  product  of  zl  and  z2 
*/ 

^define  cmul(z,  zl,  z2)\ 
double  tempi,  temp2;\ 

\ 

tempi  =  Re(zl)  *  Im(z2);\ 
temp2  =  Im(zl)  *  Re(z2);\ 

Re(z)  =  (Re(zl)  +  Im(zl))  *  (Re(z2)  -  Im(z2))  +  tempi  -  temp2;\ 
Im(z)  =  tempi  +  temp2;\ 

} 

/* 

*  z  is  the  principal  square  root  of  zl 
*/ 

Mefine  csqrtfz,  zl)\ 

{\ 

double  a,  b,  r;\ 

\ 

a  =  Re(zl);\ 
b  =  Im(zl);\ 
r  =  sqrt(a  *  a  +  b  *  b);\ 

/*\ 

*  check  for  truncation  error\ 
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*/\ 

if  ((Re(z)  =  (r  +  a)  /  2)  >  0)\ 
Re(z)  =  sqrt(Re(z));\ 
else\ 

Re(z)  =  0;\ 

if  ((Im(z)  =  (r  -  a)  /  2)  >  0)\ 
Im(z)  =  sqrt(Im(z));\ 
else\ 

Im(z)  =  0;\ 

/*\ 

*  select  the  principal  branch\ 

*/\ 

if  (b  <  0)\ 

Im(z)  =  -Im(z);\ 

} 


/* 

*  z  is  the  complex  difference  of  zl  and  z2 
*/ 

^define  csub(z,  zl,  z2)\ 

Re(z)  =  Re(zl)  -  Re(z2);\ 

Im(z)  =  Im(zl)  -  Im(z2) 


#endif 


FFP.H 


#ifndef_COMPLEX_ 

#defme  _COMPLEX_ 
typedef  struct 
{ 

double  real,  imaginary; 

}  complex; 

#endif 

struct  geometry 

{ 

double  delta_theta,  range,  theta_l,  theta_2,  Zr,  Zs; 

}; 


struct  ground 

{ 

double  depth,  omega_l,  omega_2,  pn_l,  pn_2, 
sf _ 1 ,  sf_2,  sigma_l,  sigma_2; 

}; 


struct  interface 

{ 

double  c,  kelvin,  pressure,  relative_humidity, 
wind  speed,  wind  theta,  z; 

}; 


struct  gp 

{ 

complex  impedance_l,  impedance_n,  k_l,  k_n; 
int  index  detector,  index  source; 

}; 

struct  layer 

{ 

complex  ki2,  kz,  yO; 
double  c,  mu,  rho,  thickness; 

}; 


Appendix  B 
Example  Cases 


Example  Case  #1 


Input  Files 


Debug  File 

133  (Sound  Speed  Prfile,  Source  Centered,  Wavenumber  Spectrum) 


Weather  File 


0.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

25.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

50.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

75.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

100.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

125.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

150.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

175.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

200.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

225.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

250.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

275.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

300.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

325.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

350.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

375.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

400.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

425.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

450.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

475.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

500.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

525.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

550.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

575.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

600.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

625.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

650.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

675.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

700.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

725.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

750.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

775.000 

294.8279 

0.00 

1.00 

0.00 

0.00 

800.000 

294.8279 

0.00 

1.00 

0.00 

0.00 
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825.000 

850.000 

875.000 

900.000 

925.000 

950.000 

975.000 

1000.000 

1025.000 

1050.000 

1075.000 

1100.000 

1125.000 

1150.000 

1175.000 

1200.000 

1225.000 

1250.000 

1275.000 

1300.000 

1325.000 

1350.000 

1375.000 

1400.000 

1425.000 

1450.000 

1475.000 

1500.000 

Geometry  File 
5.0 

Frequency  File 
100.0 

Ground  File 
1 

366.0 


294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 

294.8279 


1.0  10000.0 


0.0 


0.27 


0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.0 

0.0 


0.5  0.5 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

1.0 
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Output  Files 


Realtime  Output: 


Input  name  of  debug  file:  debug2 

Input  name  of  weather  file:  easel. wea 

Input  name  of  geometry  file:  geometry 2 

Input  name  of  frequency  file:  frequency2 

Input  name  of  ground  file:  ground2 


SINGLE  FREQUENCIES 


RECEIVER  HEIGHT  (Zr) 

=  1.0  METERS 

SOURCE  HEIGHT  (Zs) 

RANGE 

=  10000.0  METERS  BEARING  OF  SOURCE 

FREQUENCY 

=  100  Hz 

KELVIN 

RELATIVE 

ATMOSPHERIC 

INTERFACE 

HEIGHT 

TEMPERATURE 

HUMIDITY 

PRESSURE 

61 

1500.00 

294.83 

0.0 

1.000 

60 

1475.00 

294.83 

0.0 

1.000 

59 

1450.00 

294.83 

0.0 

1.000 

58 

1425.00 

294.83 

0.0 

1.000 

57 

1400.00 

294.83 

0.0 

1.000 

56 

1375.00 

294.83 

0.0 

1.000 

55 

1350.00 

294.83 

0.0 

1.000 

54 

1325.00 

294.83 

0.0 

1.000 

53 

1300.00 

294.83 

0.0 

1.000 

52 

1275.00 

294.83 

0.0 

1.000 

51 

1250.00 

294.83 

0.0 

1.000 

50 

1225.00 

294.83 

0.0 

1.000 

49 

1200.00 

294.83 

0.0 

1.000 

48 

1175.00 

294.83 

0.0 

1.000 

47 

1150.00 

294.83 

0.0 

1.000 

46 

1125.00 

294.83 

0.0 

1.000 

45 

1100.00 

294.83 

0.0 

1.000 

44 

1075.00 

294.83 

0.0 

1.000 

43 

1050.00 

294.83 

0.0 

1.000 

42 

1025.00 

294.83 

0.0 

1.000 

41 

1000.00 

294.83 

0.0 

1.000 

40 

975.00 

294.83 

0.0 

1.000 

39 

950.00 

294.83 

0.0 

1.000 

38 

925.00 

294.83 

0.0 

1.000 

37 

900.00 

294.83 

0.0 

1.000 

=  5.0  METERS 

=  0.000000  DEGREES 


WIND  WIND 

SPEED  DIRECTION 
0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  o.ooo 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 
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36 

875.00 

294.83 

35 

850.00 

294.83 

34 

825.00 

294.83 

33 

800.00 

294.83 

32 

775.00 

294.83 

31 

750.00 

294.83 

30 

725.00 

294.83 

29 

700.00 

294.83 

28 

675.00 

294.83 

27 

650.00 

294.83 

26 

625.00 

294.83 

25 

600.00 

294.83 

24 

575.00 

294.83 

23 

550.00 

294.83 

22 

525.00 

294.83 

21 

500.00 

294.83 

20 

475.00 

294.83 

19 

450.00 

294.83 

18 

425.00 

294.83 

17 

400.00 

294.83 

16 

375.00 

294.83 

15 

350.00 

294.83 

14 

325.00 

294.83 

13 

300.00 

294.83 

12 

275.00 

294.83 

11 

250.00 

294.83 

10 

225.00 

294.83 

9 

200.00 

294.83 

8 

175.00 

294.83 

7 

150.00 

294.83 

6 

125.00 

294.83 

5 

100.00 

294.83 

4 

75.00 

294.83 

3 

50.00 

294.83 

2 

25.00 

294.83 

1 

0.00 

294.83 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.0 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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GROUND  LAYER  SIGMA  POROSITY 

1  366.00  0.27 


Position  Angle  =  0.000000 
10000.00  -124.9824 


Spectrum 


Case  #1 


Example  Case  #2 


Input  Files 


Debug  File 
133 


Weather  File 
61 

0.000 

25.000 

50.000 

75.000 

100.000 

125.000 

150.000 

175.000 

200.000 

225.000 

250.000 

275.000 

300.000 

325.000 

350.000 

375.000 

400.000 

425.000 

450.000 

475.000 

500.000 

525.000 

550.000 

575.000 

600.000 

625.000 

650.000 

675.000 

700.000 

725.000 

750.000 

775.000 

800.000 


(Sound  Speed  Prfile,  Source  Centered,  Wavenumber  Spectrum) 


294.8279 

0.00 

1.00 

0.00 

299.1287 

0.00 

1.00 

0.00 

303.4607 

0.00 

1.00 

0.00 

307.8239 

0.00 

1.00 

0.00 

312.2182 

0.00 

1.00 

0.00 

316.6436 

0.00 

1.00 

0.00 

321.1002 

0.00 

1.00 

0.00 

325.5879 

0.00 

1.00 

0.00 

330.1067 

0.00 

1.00 

0.00 

334.6567 

0.00 

1.00 

0.00 

339.2379 

0.00 

1.00 

0.00 

343.8502 

0.00 

1.00 

0.00 

348.4936 

0.00 

1.00 

0.00 

353.1682 

0.00 

1.00 

0.00 

357.8739 

0.00 

1.00 

0.00 

362.6108 

0.00 

1.00 

0.00 

367.3788 

0.00 

1.00 

0.00 

372.1779 

0.00 

1.00 

0.00 

377.0082 

0.00 

1.00 

0.00 

381.8697 

0.00 

1.00 

0.00 

386.7622 

0.00 

1.00 

0.00 

391.6860 

0.00 

1.00 

0.00 

396.6408 

0.00 

1.00 

0.00 

401.6268 

0.00 

1.00 

0.00 

406.6440 

0.00 

1.00 

0.00 

411.6923 

0.00 

1.00 

0.00 

416.7717 

0.00 

1.00 

0.00 

421.8823 

0.00 

1.00 

0.00 

427.0240 

0.00 

1.00 

0.00 

432.1969 

0.00 

1.00 

0.00 

437.4009 

0.00 

1.00 

0.00 

442.6361 

0.00 

1.00 

0.00 

447.9024 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 
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825.000 

453.1998 

0.00 

1.00 

0.00 

0.00 

850.000 

458.5284 

0.00 

1.00 

0.00 

0.00 

875.000 

463.8881 

0.00 

1.00 

0.00 

0.00 

900.000 

469.2790 

0.00 

1.00 

0.00 

0.00 

925.000 

474.7010 

0.00 

1.00 

0.00 

0.00 

950.000 

480.1541 

0.00 

1.00 

0.00 

0.00 

975.000 

485.6384 

0.00 

1.00 

0.00 

0.00 

1000.000 

491.1539 

0.00 

1.00 

0.00 

0.00 

1025.000 

496.7005 

0.00 

1.00 

0.00 

0.00 

1050.000 

502.2782 

0.00 

1.00 

0.00 

0.00 

1075.000 

507.8871 

0.00 

1.00 

0.00 

0.00 

1100.000 

513.5271 

0.00 

1.00 

0.00 

0.00 

1125.000 

519.1982 

0.00 

1.00 

0.00 

0.00 

1150.000 

524.9005 

0.00 

1.00 

0.00 

0.00 

1175.000 

530.6340 

0.00 

1.00 

0.00 

0.00 

1200.000 

536.3985 

0.00 

1.00 

0.00 

0.00 

1225.000 

542.1943 

0.00 

1.00 

0.00 

0.00 

1250.000 

548.0211 

0.00 

1.00 

0.00 

0.00 

1275.000 

553.8792 

0.00 

1.00 

0.00 

0.00 

1300.000 

559.7683 

0.00 

1.00 

0.00 

0.00 

1325.000 

565.6886 

0.00 

1.00 

0.00 

0.00 

1350.000 

571.6401 

0.00 

1.00 

0.00 

0.00 

1375.000 

577.6227 

0.00 

1.00 

0.00 

0.00 

1400.000 

583.6364 

0.00 

1.00 

0.00 

0.00 

1425.000 

589.6813 

0.00 

1.00 

0.00 

0.00 

1450.000 

595.7573 

0.00 

1.00 

0.00 

0.00 

1475.000 

601.8644 

0.00 

1.00 

0.00 

0.00 

1500.000 

608.0027 

0.00 

1.00 

0.00 

0.00 

Geometry  File 

5.0 

1.0 

10000.0 

0.0 

0.0 

1.0 

Frequency  File 

100.0 

Ground  File 

1 

366.0 

0.0 

0.0 

0.27 

0.5 

0.5 
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Output  Files 


Realtime  Output: 


Input  name  of  debug  file:  debug2 

Input  name  of  weather  file:  case2.wea 

Input  name  of  geometry  file:  geometry2 

Input  name  of  frequency  file:  frequency2 

Input  name  of  ground  fde:  ground2 


RECEIVER  HEIGHT  (Zr)  =  1.0  METERS 

RANGE  =  10000.0  METERS 

FREQUENCY  =  100  Hz 


SINGLE  FREQUENCIES 


SOURCE  HEIGHT  (Zs)  =  5.0  METERS 

BEARING  OF  SOURCE  =  0.000000  DEGREES 


KELVIN 

RELATIVE 

ATMOSPHERIC 

WIND 

WIND 

INTERFACE 

HEIGHT 

TEMPERATURE 

HUMIDITY 

PRESSURE 

SPEED 

DIRECTION 

61 

1500.00 

608.00 

0.0 

1.000 

0.000 

0.000 

60 

1475.00 

601.86 

0.0 

1.000 

0.000 

0.000 

59 

1450.00 

595.76 

0.0 

1.000 

0.000 

0.000 

58 

1425.00 

589.68 

0.0 

1.000 

0.000 

0.000 

57 

1400.00 

583.64 

0.0 

1.000 

0.000 

0.000 

56 

1375.00 

577.62 

0.0 

1.000 

0.000 

0.000 

55 

1350.00 

571.64 

0.0 

1.000 

0.000 

0.000 

54 

1325.00 

565.69 

0.0 

1.000 

0.000 

0.000 

53 

1300.00 

559.77 

0.0 

1.000 

0.000 

0.000 

52 

1275.00 

553.88 

0.0 

1.000 

0.000 

0.000 

51 

1250.00 

548.02 

0.0 

1.000 

0.000 

0.000 

50 

1225.00 

542.19 

0.0 

1.000 

0.000 

0.000 

49 

1200.00 

536.40 

0.0 

1.000 

0.000 

0.000 

48 

1175.00 

530.63 

0.0 

1.000 

0.000 

0.000 

47 

1150.00 

524.90 

0.0 

1.000 

0.000 

0.000 

46 

1125.00 

519.20 

0.0 

1.000 

0.000 

0.000 

45 

1100.00 

513.53 

0.0 

1.000 

0.000 

0.000 

44 

1075.00 

507.89 

0.0 

1.000 

0.000 

0.000 

43 

1050.00 

502.28 

0.0 

1.000 

0.000 

0.000 

42 

1025.00 

496.70 

0.0 

1.000 

0.000 

0.000 

41 

1000.00 

491.15 

0.0 

1.000 

0.000 

0.000 

40 

975.00 

485.64 

0.0 

1.000 

0.000 

0.000 

39 

950.00 

480.15 

0.0 

1.000 

0.000 

0.000 

38 

925.00 

474.70 

0.0 

1.000 

0.000 

0.000 

37 

900.00 

469.28 

0.0 

1.000 

0.000 

0.000 
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36 

35 

34 

33 

32 

31 

30 

29 

28 

27 

26 

25 

24 

23 

22 

21 

20 

19 

18 

17 

16 

15 

14 

13 

12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 


875.00 

463.89 

0.0 

850.00 

458.53 

0.0 

825.00 

453.20 

0.0 

800.00 

447.90 

0.0 

775.00 

442.64 

0.0 

750.00 

437.40 

0.0 

725.00 

432.20 

0.0 

700.00 

427.02 

0.0 

675.00 

421.88 

0.0 

650.00 

416.77 

0.0 

625.00 

411.69 

0.0 

600.00 

406.64 

0.0 

575.00 

401.63 

0.0 

550.00 

396.64 

0.0 

525.00 

391.69 

0.0 

500.00 

386.76 

0.0 

475.00 

381.87 

0.0 

450.00 

377.01 

0.0 

425.00 

372.18 

0.0 

400.00 

367.38 

0.0 

375.00 

362.61 

0.0 

350.00 

357.87 

0.0 

325.00 

353.17 

0.0 

300.00 

348.49 

0.0 

275.00 

343.85 

0.0 

250.00 

339.24 

0.0 

225.00 

334.66 

0.0 

200.00 

330.11 

0.0 

175.00 

325.59 

0.0 

150.00 

321.10 

0.0 

125.00 

316.64 

0.0 

100.00 

312.22 

0.0 

75.00 

307.82 

0.0 

50.00 

303.46 

0.0 

25.00 

299.13 

0.0 

0.00 

294.83 

0.0 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 

1.000 

0.000 

0.000 
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GROUND  LAYER  SIGMA  POROSITY 

1  366.00  0.27 


Position  Angle  =  0.000000 
10000.00  -96.3082 


Sound  Speed  Profile: 


Case  #2 


94 


Magnitude  of  Spectrum  |  Relative  Sound  Pressure  Level  (dB) 


Relative  Sound  Pressure  Level  versus  Range: 


Case  #2 


enumber  Spectrum: 


Case  #2 


95 


Example  Case  #3 


Input  Files 


Debug  File 
133 


Weather  File 
61 

0.000 

25.000 

50.000 

75.000 

100.000 

125.000 

150.000 

175.000 

200.000 

225.000 

250.000 

275.000 

300.000 

325.000 

350.000 

375.000 

400.000 

425.000 

450.000 

475.000 

500.000 

525.000 

550.000 

575.000 

600.000 

625.000 

650.000 

675.000 

700.000 

725.000 

750.000 

775.000 

800.000 


(Sound  Speed  Prfde,  Source  Centered,  Wavenumber  Spectrum) 


294.8279 

0.00 

1.00 

0.00 

290.5581 

0.00 

1.00 

0.00 

286.3196 

0.00 

1.00 

0.00 

282.1121 

0.00 

1.00 

0.00 

277.9358 

0.00 

1.00 

0.00 

273.7907 

0.00 

1.00 

0.00 

269.6767 

0.00 

1.00 

0.00 

265.5938 

0.00 

1.00 

0.00 

261.5421 

0.00 

1.00 

0.00 

257.5216 

0.00 

1.00 

0.00 

253.5321 

0.00 

1.00 

0.00 

249.5738 

0.00 

1.00 

0.00 

245.6467 

0.00 

1.00 

0.00 

241.7507 

0.00 

1.00 

0.00 

237.8858 

0.00 

1.00 

0.00 

234.0521 

0.00 

1.00 

0.00 

230.2496 

0.00 

1.00 

0.00 

226.4781 

0.00 

1.00 

0.00 

222.7378 

0.00 

1.00 

0.00 

219.0287 

0.00 

1.00 

0.00 

215.3507 

0.00 

1.00 

0.00 

211.7038 

0.00 

1.00 

0.00 

208.0881 

0.00 

1.00 

0.00 

204.5036 

0.00 

1.00 

0.00 

200.9501 

0.00 

1.00 

0.00 

197.4279 

0.00 

1.00 

0.00 

193.9367 

0.00 

1.00 

0.00 

190.4767 

0.00 

1.00 

0.00 

187.0479 

0.00 

1.00 

0.00 

183.6502 

0.00 

1.00 

0.00 

180.2836 

0.00 

1.00 

0.00 

176.9482 

0.00 

1.00 

0.00 

173.6439 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 
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825.000 

170.3707 

0.00 

1.00 

0.00 

0.00 

850.000 

167.1287 

0.00 

1.00 

0.00 

0.00 

875.000 

163.9179 

0.00 

1.00 

0.00 

0.00 

900.000 

160.7382 

0.00 

1.00 

0.00 

0.00 

925.000 

157.5896 

0.00 

1.00 

0.00 

0.00 

950.000 

154.4722 

0.00 

1.00 

0.00 

0.00 

975.000 

151.3859 

0.00 

1.00 

0.00 

0.00 

1000.000 

148.3308 

0.00 

1.00 

0.00 

0.00 

1025.000 

145.3068 

0.00 

1.00 

0.00 

0.00 

1050.000 

142.3139 

0.00 

1.00 

0.00 

0.00 

1075.000 

139.3522 

0.00 

1.00 

0.00 

0.00 

1100.000 

136.4217 

0.00 

1.00 

0.00 

0.00 

1125.000 

133.5222 

0.00 

1.00 

0.00 

0.00 

1150.000 

130.6540 

0.00 

1.00 

0.00 

0.00 

1175.000 

127.8168 

0.00 

1.00 

0.00 

0.00 

1200.000 

125.0108 

0.00 

1.00 

0.00 

0.00 

1225.000 

122.2360 

0.00 

1.00 

0.00 

0.00 

1250.000 

119.4923 

0.00 

1.00 

0.00 

0.00 

1275.000 

116.7797 

0.00 

1.00 

0.00 

0.00 

1300.000 

114.0983 

0.00 

1.00 

0.00 

0.00 

1325.000 

111.4480 

0.00 

1.00 

0.00 

0.00 

1350.000 

108.8289 

0.00 

1.00 

0.00 

0.00 

1375.000 

106.2409 

0.00 

1.00 

0.00 

0.00 

1400.000 

103.6841 

0.00 

1.00 

0.00 

0.00 

1425.000 

101.1583 

0.00 

1.00 

0.00 

0.00 

1450.000 

98.6638 

0.00 

1.00 

0.00 

0.00 

1475.000 

96.2004 

0.00 

1.00 

0.00 

0.00 

1500.000 

93.7681 

0.00 

1.00 

0.00 

0.00 

Geometry  File 

5.0 

1.0 

10000.0 

0.0 

0.0 

1.0 

Frequency  File 

100.0 

Ground  File 

1 

366.0 

0.0 

0.0 

0.27 

0.5 

0.5 
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Output  Files 


Realtime  Output: 

Input  name  of  debug  file: 

debug2 

Input  name  of  weather  file: 

case3.wea 

Input  name  of  geometry  file: 

geometry2 

Input  name  of  frequency  file: 

frequency2 

Input  name  of  ground  file: 

ground2 

RECEIVER  HEIGHT  (Zr)  =  1.0  METERS 
RANGE  =  10000.0  METERS 

FREQUENCY  =  100  Hz 


SINGLE  FREQUENCIES 


SOURCE  HEIGHT  (Zs)  =  5.0  METERS 

BEARING  OF  SOURCE  =  0.000000  DEGREES 


KELVIN 

RELATIVE 

ATMOSPHERIC 

WIND 

WIND 

INTERFACE 

HEIGHT 

TEMPERATURE 

HUMIDITY 

PRESSURE 

SPEED 

DIRECTION 

61 

1500.00 

93.77 

0.0 

1.000 

0.000 

0.000 

60 

1475.00 

96.20 

0.0 

1.000 

0.000 

0.000 

59 

1450.00 

98.66 

0.0 

1.000 

0.000 

0.000 

58 

1425.00 

101.16 

0.0 

1.000 

0.000 

0.000 

57 

1400.00 

103.68 

0.0 

1.000 

0.000 

0.000 

56 

1375.00 

106.24 

0.0 

1.000 

0.000 

0.000 

55 

1350.00 

108.83 

0.0 

1.000 

0.000 

0.000 

54 

1325.00 

111.45 

0.0 

1.000 

0.000 

0.000 

53 

1300.00 

114.10 

0.0 

1.000 

0.000 

0.000 

52 

1275.00 

116.78 

0.0 

1.000 

0.000 

0.000 

51 

1250.00 

119.49 

0.0 

1.000 

0.000 

0.000 

50 

1225.00 

122.24 

0.0 

1.000 

0.000 

0.000 

49 

1200.00 

125.01 

0.0 

1.000 

0.000 

0.000 

48 

1175.00 

127.82 

0.0 

1.000 

0.000 

0.000 

47 

1150.00 

130.65 

0.0 

1.000 

0.000 

0.000 

46 

1125.00 

133.52 

0.0 

1.000 

0.000 

0.000 

45 

1100.00 

136.42 

0.0 

1.000 

0.000 

0.000 

44 

1075.00 

139.35 

0.0 

1.000 

0.000 

0.000 

43 

1050.00 

142.31 

0.0 

1.000 

0.000 

0.000 

42 

1025.00 

145.31 

0.0 

1.000 

0.000 

0.000 

41 

1000.00 

148.33 

0.0 

1.000 

0.000 

0.000 

40 

975.00 

151.39 

0.0 

1.000 

0.000 

0.000 

39 

950.00 

154.47 

0.0 

1.000 

0.000 

0.000 

38 

925.00 

157.59 

0.0 

1.000 

0.000 

0.000 

37 

900.00 

160.74 

0.0 

1.000 

0.000 

0.000 

98 


36 

875.00 

163.92 

0.0 

1.000 

0.000 

35 

850.00 

167.13 

0.0 

1.000 

0.000 

34 

825.00 

170.37 

0.0 

1.000 

0.000 

33 

800.00 

173.64 

0.0 

1.000 

0.000 

32 

775.00 

176.95 

0.0 

1.000 

0.000 

31 

750.00 

180.28 

0.0 

1.000 

0.000 

30 

725.00 

183.65 

0.0 

1.000 

0.000 

29 

700.00 

187.05 

0.0 

1.000 

0.000 

28 

675.00 

190.48 

0.0 

1.000 

0.000 

27 

650.00 

193.94 

0.0 

1.000 

0.000 

26 

625.00 

197.43 

0.0 

1.000 

0.000 

25 

600.00 

200.95 

0.0 

1.000 

0.000 

24 

575.00 

204.50 

0.0 

1.000 

0.000 

23 

550.00 

208.09 

0.0 

1.000 

0.000 

22 

525.00 

211.70 

0.0 

1.000 

0.000 

21 

500.00 

215.35 

0.0 

1.000 

0.000 

20 

475.00 

219.03 

0.0 

1.000 

0.000 

19 

450.00 

222.74 

0.0 

1.000 

0.000 

18 

425.00 

226.48 

0.0 

1.000 

0.000 

17 

400.00 

230.25 

0.0 

1.000 

0.000 

16 

375.00 

234.05 

0.0 

1.000 

0.000 

15 

350.00 

237.89 

0.0 

1.000 

0.000 

14 

325.00 

241.75 

0.0 

1.000 

0.000 

13 

300.00 

245.65 

0.0 

1.000 

0.000 

12 

275.00 

249.57 

0.0 

1.000 

0.000 

11 

250.00 

253.53 

0.0 

1.000 

0.000 

10 

225.00 

257.52 

0.0 

1.000 

0.000 

9 

200.00 

261.54 

0.0 

1.000 

0.000 

8 

175.00 

265.59 

0.0 

1.000 

0.000 

7 

150.00 

269.68 

0.0 

1.000 

0.000 

6 

125.00 

273.79 

0.0 

1.000 

0.000 

5 

100.00 

277.94 

0.0 

1.000 

0.000 

4 

75.00 

282.11 

0.0 

1.000 

0.000 

3 

50.00 

286.32 

0.0 

1.000 

0.000 

2 

25.00 

290.56 

0.0 

1.000 

0.000 

1 

0.00 

294.83 

0.0 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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GROUND  LAYER  SIGMA  POROSITY 

1  366.00  0.27 

Position  Angle  =  0.000000 
10000.00  -143.1275 


Example  Case  #4 


Debug  File 
133 

Weather  File 
65 

0.000 

5.000 

10.000 

15.000 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 

55.000 

60.000 

65.000 

70.000 

75.000 

80.000 

85.000 

90.000 

95.000 

100.000 

110.000 

120.000 

130.000 

140.000 

150.000 

160.000 

170.000 

180.000 

190.000 

200.000 

210.000 

220.000 


Input  Files 

(Sound  Speed  Prfile,  Source  Centered,  Wavenumber  Spectrum) 


294.8279 

0.00 

1.00 

0.00 

0.00 

295.6855 

0.00 

1.00 

0.00 

0.00 

296.5445 

0.00 

1.00 

0.00 

0.00 

297.4046 

0.00 

1.00 

0.00 

0.00 

298.2661 

0.00 

1.00 

0.00 

0.00 

299.1287 

0.00 

1.00 

0.00 

0.00 

299.9926 

0.00 

1.00 

0.00 

0.00 

300.8578 

0.00 

1.00 

0.00 

0.00 

301.7242 

0.00 

1.00 

0.00 

0.00 

302.5918 

0.00 

1.00 

0.00 

0.00 

303.4607 

0.00 

1.00 

0.00 

0.00 

304.3309 

0.00 

1.00 

0.00 

0.00 

305.2022 

0.00 

1.00 

0.00 

0.00 

306.0749 

0.00 

1.00 

0.00 

0.00 

306.9487 

0.00 

1.00 

0.00 

0.00 

307.8239 

0.00 

1.00 

0.00 

0.00 

308.7002 

0.00 

1.00 

0.00 

0.00 

309.5778 

0.00 

1.00 

0.00 

0.00 

310.4567 

0.00 

1.00 

0.00 

0.00 

311.3368 

0.00 

1.00 

0.00 

0.00 

312.2182 

0.00 

1.00 

0.00 

0.00 

310.4567 

0.00 

1.00 

0.00 

0.00 

308.7002 

0.00 

1.00 

0.00 

0.00 

306.9487 

0.00 

1.00 

0.00 

0.00 

305.2022 

0.00 

1.00 

0.00 

0.00 

303.4607 

0.00 

1.00 

0.00 

0.00 

301.7242 

0.00 

1.00 

0.00 

0.00 

299.9926 

0.00 

1.00 

0.00 

0.00 

298.2661 

0.00 

1.00 

0.00 

0.00 

296.5445 

0.00 

1.00 

0.00 

0.00 

294.8279 

0.00 

1.00 

0.00 

0.00 

293.1162 

0.00 

1.00 

0.00 

0.00 

291.4096 

0.00 

1.00 

0.00 

0.00 
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230.000 

289.7079 

0.00 

1.00 

0.00 

0.00 

240.000 

288.0113 

0.00 

1.00 

0.00 

0.00 

250.000 

286.3196 

0.00 

1.00 

0.00 

0.00 

260.000 

284.6329 

0.00 

1.00 

0.00 

0.00 

270.000 

282.9511 

0.00 

1.00 

0.00 

0.00 

280.000 

281.2744 

0.00 

1.00 

0.00 

0.00 

290.000 

279.6026 

0.00 

1.00 

0.00 

0.00 

300.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

350.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

400.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

450.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

500.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

550.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

600.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

650.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

700.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

750.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

800.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

850.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

900.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

950.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1000.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1050.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1100.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1150.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1200.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1250.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1300.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1350.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1400.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1450.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

1500.000 

277.9358 

0.00 

1.00 

0.00 

0.00 

Geometry  File 

5.0 

1.0 

10000.0 

0.0 

0.0 

1.0 

Frequency  File 

100.0 

0.0 

0.0 

Ground  File 

1 

366.0 

0.27 

0.5 

0.5 
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Output  Files 


Realtime  Output: 


Input  name  of  debug  file:  debug2 

Input  name  of  weather  file:  case4.wea 

Input  name  of  geometry  file:  geometry2 

Input  name  of  frequency  file:  frequency2 

Input  name  of  ground  file:  ground2 


RECEIVER  HEIGHT  (Zr) 

RANGE 

FREQUENCY 


rERFACE 

HEIGHT 

65 

1500.00 

64 

1450.00 

63 

1400.00 

62 

1350.00 

61 

1300.00 

60 

1250.00 

59 

1200.00 

58 

1150.00 

57 

1100.00 

56 

1050.00 

55 

1000.00 

54 

950.00 

53 

900.00 

52 

850.00 

51 

800.00 

50 

750.00 

49 

700.00 

48 

650.00 

47 

600.00 

46 

550.00 

45 

500.00 

44 

450.00 

43 

400.00 

42 

350.00 

41 

300.00 

40 

290.00 

39 

280.00 

38 

270.00 

37 

260.00 

=  1.0  METERS 
=  10000.0  METERS 
=  100  Hz 


KELVIN 

TEMPERATURE 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 

277.94 
279.60 
281.27 

282.95 
284.63 


SINGLE  FREQUENCIES 


SOURCE  HEIGHT  (Zs) 
BEARING  OF  SOURCE 


RELATIVE 

HUMIDITY 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


ATMOSPHERIC 

PRESSURE 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


=  5.0  METERS 

=  0.000000  DEGREES 


WIND 

WIND 

SPEED 

DIRECTIC 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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36 

250.00 

286.32 

0.0 

1.000 

0.000 

0.000 

35 

240.00 

288.01 

0.0 

1.000 

0.000 

0.000 

34 

230.00 

289.71 

0.0 

1.000 

0.000 

0.000 

33 

220.00 

291.41 

0.0 

1.000 

0.000 

0.000 

32 

210.00 

293.12 

0.0 

1.000 

0.000 

0.000 

31 

200.00 

294.83 

0.0 

1.000 

0.000 

0.000 

30 

190.00 

296.54 

0.0 

1.000 

0.000 

0.000 

29 

180.00 

298.27 

0.0 

1.000 

0.000 

0.000 

28 

170.00 

299.99 

0.0 

1.000 

0.000 

0.000 

27 

160.00 

301.72 

0.0 

1.000 

0.000 

0.000 

26 

150.00 

303.46 

0.0 

1.000 

0.000 

0.000 

25 

140.00 

305.20 

0.0 

1.000 

0.000 

0.000 

24 

130.00 

306.95 

0.0 

1.000 

0.000 

0.000 

23 

120.00 

308.70 

0.0 

1.000 

0.000 

0.000 

22 

110.00 

310.46 

0.0 

1.000 

0.000 

0.000 

21 

100.00 

312.22 

0.0 

1.000 

0.000 

0.000 

20 

95.00 

311.34 

0.0 

1.000 

0.000 

0.000 

19 

90.00 

310.46 

0.0 

1.000 

0.000 

0.000 

18 

85.00 

309.58 

0.0 

1.000 

0.000 

0.000 

17 

80.00 

308.70 

0.0 

1.000 

0.000 

0.000 

16 

75.00 

307.82 

0.0 

1.000 

0.000 

0.000 

15 

70.00 

306.95 

0.0 

1.000 

0.000 

0.000 

14 

65.00 

306.07 

0.0 

1.000 

0.000 

0.000 

13 

60.00 

305.20 

0.0 

1.000 

0.000 

0.000 

12 

55.00 

304.33 

0.0 

1.000 

0.000 

0.000 

11 

50.00 

303.46 

0.0 

1.000 

0.000 

0.000 

10 

45.00 

302.59 

0.0 

1.000 

0.000 

0.000 

9 

40.00 

301.72 

0.0 

1.000 

0.000 

0.000 

8 

35.00 

300.86 

0.0 

1.000 

0.000 

0.000 

7 

30.00 

299.99 

0.0 

1.000 

0.000 

0.000 

6 

25.00 

299.13 

0.0 

1.000 

0.000 

0.000 

5 

20.00 

298.27 

0.0 

1.000 

0.000 

0.000 

4 

15.00 

297.40 

0.0 

1.000 

0.000 

0.000 

3 

10.00 

296.54 

0.0 

1.000 

0.000 

0.000 

2 

5.00 

295.69 

0.0 

1.000 

0.000 

0.000 

1 

0.00 

294.83 

0.0 

1.000 

0.000 

0.000 
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GROUND  LAYER  SIGMA  POROSITY 

1  366.00  0.27 

Position  Angle  =  0.000000 
10000.00  -114.5068 


Magnitude  of  Spectrum  f  Relative  Sound  Pressure  Level  (dB) 


Relative  Sound  Pressure  Level  versus  Range: 


Case  #4 


number  Spectrum: 


Case  #4 
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Appendix  C 
Debug  Features 


As  mentioned  in  section  3.2,  SCAFFIP  allows  for  input  of  a  debug  code  to 
allow  for  modification  of  some  of  the  default  parameters  and  provide  some 
additional  output  information.  This  allows  for  a  standardization  of  the  input 
and  output  format,  yet  allow  for  versatility.  Currently,  the  debug  files  consist 
of  up  to  six  files:  two  input  files  and  four  output  files. 

Input  Files 

Code:  8 

File  Name:  extra. los 
Format: 
extra 

where 

extra  -  value  for  extra  loss  typically  between  10‘6  to  10~2  (default 

10-4) 

Code:  16 

File  Name:  trans.frm 
Format: 

points  Npan 
where 

points  -  number  of  points  to  use  in  calculations.  Must  be 
power  of  2.  (default  1024) 

Npan  -  set  the  number  of  panels  to  use  in  calculation. 

Npan  <  0  Calculate  number  of  panels  (default  -2) 
Npan  >  0  Set  number  of  panels 
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Output  Files 


Code:  1 

File  Name: 

sound,  o 

Format: 

Z(l) 

C(l) 

Z(2) 

C(2) 

Z(Nint; 

)  C(Nint) 

where 

Nint 

number  of  interfaces. 

Z(i) 

height  of  the  ith  interface  in  meters. 

C(i) 

sound  speed  of  the  ith  interface  in  m/s. 

Code:  32 

File  Name: 

npan.e 

Format: 

Kmax 

Nyq  Npan 

Kwidth 

DelR  Delk 

where 

Kmax 

maximum  wave  number  in  summation 

DelR 

range  increment 

Delk 

wavenumber  increment 

Npan 

number  of  panels 

Kwidth 

wavenumber  band  width  per  panel 

Nyq 

nyquist  criteria 

Code:  64 

File  Name: 

imped,  e 

Format: 

Zreal 

Zimg 

Kreal 

Kimg 

where 


Zreal,  Zimg  -  real  and  imaginary  part  of  the  ground  impedance 
Kreal,  Kimg  -  real  and  imaginary  part  of  the  wavenumber 

Code:  128 

File  Name:  wavnum.e 
Format: 

wavenum(l)  mag(l) 

wavenum(2)  mag(2) 


wavenum(points)  mag(points) 


where 

points 

wavnum(i) 

mag(i) 


number  of  points  in  transform 
wavenumber  of  ith  point 
magnitude  of  the  kernal  of  the  ith  point 


Appendix  D 

User  Developed  Interface 
Requirements  for  SCAFFIP 


This  appendix  discusses  the  possibility  of  the  user  wishing  to  develop  a 
personalized  front  end  user  interface  for  SC  AFFIP .  As  mentioned  in  section  1 , 
deviations  from  the  standard  user  interface  supported  by  ARL/BED  will  not 
automatically  be  provided  debugging  assistance.  Future  upgrades  will  conform 
to  the  standard  interface  and  will  not  be  integrated  into  a  user  developed 
interface.  This  policy  is  in  effect  because  of  funding  and  personnel  shortages. 
Limited  assistance  can  be  provided  for  input  data  and  problems  related  to  the 
propagation  models .  Appendix  A  contains  the  complete  source  code  listing  for 

the  model. 

The  front-end  user  interface  for  SCAFFIP  is  configured  so  the  five  input  files 
are  read.  There  are  some  additional  input  files  that  may  be  read;  however, 
these  files  deal  with  the  debug  number  read  in.  The  five  primary  input  files 
are  discussed  here.  Debug  input  files  are  discussed  in  appendix  C.  The  front- 
end  user  interface  files  for  SCAFFIP  are  main.c,  input. c,  header. c,  parser. c, 
and  sound. c.  The  file  main.c  is  the  main  driver  program  that  calls  all  the  other 
parts  of  the  program.  The  file  input. c  is  the  function  that  reads  the  five  input 
files  and  calls  parser  to  decipher  the  debug  code.  The  file  header.c  provides 
a  visual  printout  of  the  input  parameters.  The  file  sound. c  calculates  the  sound 
speed  profile  from  the  meteorological  data  and  the  given  geometry.  Because 
the  sound  speed  profile  is  dependent  on  a  vector  relationship  between  the  wind 
direction  and  the  direction  of  propagation,  the  sound  speed  profile  is  calculated 
at  each  azimuth  value.  The  functions  input  and  sound  initialize  three  of  the 
structures  contained  in  the  include  file  ffp.h. 

To  develop  a  new  user  interface,  the  new  user  interface  must  initialize  these 
three  structures:  (1)  geometry,  (2)  ground,  and  (3)  interface.  The  include  file 
ffp.h  contains  two  other  structures;  however,  these  structures  are  initialized  and 
used  within  propmod  exclusively.  The  geometry  structure  contains  the 
parameters  that  describe  the  structure  of  the  geometry  of  the  problem:  source 
height,  receiver  height,  horizontal  range,  starting  scan  angle,  ending  scan 
angle,  and  angular  resolution  of  scan.  The  ground  structure  contains  all  of  the 
parameters  for  calculating  the  ground  impedance  using  Attenborough’s  Four 
Parameter  Model  (see  section  2.4).  The  interface  structure  contains  the 
meteorological  profile  data.  Although  this  data  is  used  in  the  function  sound, 
some  of  this  information  is  required  by  the  function  propmod  to  calculate  the 


molecular  absorption  coefficients  (see  section  2.3).  While  modifying  the 
standard  user  interface,  the  user  may  choose  to  input  the  values  for  the  ground 
impedance  directly  instead  of  using  the  built-in  ground  impedance  model. 
This  can  be  accomplished  by  storing  the  values  into  the  structure  gp  and 
replacing  the  variable  ground  in  the  argument  list  for  propmod  by  a  variable 
of  type  gp  that  contains  the  user  input  ground  impedance  value.  An  alternative 
method  for  storing  the  ground  impedance  is  to  use  one  of  the  unused  debug 
numbers  and  place  a  conditional  statement  to  read  the  file  containing  the 
ground  impedance  values  or  use  the  ground  impedance  model  to  calculate  the 
values.  The  second  method  is  preferred  because  it  conforms  better  to  the 
programming  philosophy  used  in  the  code  development. 

After  the  three  structures  are  initialized,  the  control  can  be  passed  off  to  the 
propagation  function  propmod.  The  argument  list  for  propmod  is 


a 

fd 

fp 

frequency 

srclevel 

backlevel 

geometry 

ground 

n_ground_layers 

ninterfaces 

dflags 


interface  structure  array 

integer  file  pointer  to  a  binary  output  file 

file  pointer  to  an  ASCII  output  file 

frequency  of  interest 

source  level 

background  level 

geometry  structure 

ground  structure 

number  of  ground  layers 

number  of  atmospheric  interfaces 

array  of  debug  parameters. 


When  the  function  propmod  is  called,  these  are  the  arguments  that  must  be 
passed.  These  are  the  areas  that  must  be  considered  when  developing  a  new 
user  interface. 


If  the  output  file  needs  to  be  modified,  the  computer  code  for  storing  the  model 
results  is  located  at  the  end  of  the  function  propmod.  The  open  statement  for 
the  output  file  is  located  in  main  because  main  controls  the  scanning  loop.  If 
the  number  of  files  opened  at  once  were  to  be  minimized,  the  open  statement 
could  be  placed  in  propmod  also;  however,  the  open  statement  would  have  to 
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be  in  setup  to  always  append  to  the  file.  This  method  would  only  allow  the 
output  file  to  be  open  when  there  is  data  to  be  written.  As  an  error  check,  a 
statement  would  have  to  be  added  to  check  for  the  presence  of  the  output  file 
at  the  beginning  of  the  program  to  prevent  appending  to  a  file  already  present. 

These  are  the  basic  areas  in  which  the  user  may  wish  to  alter  the  program. 
Other  areas  may  deal  with  the  input  and  output  of  the  parameters  associated 
with  the  debug  parameters.  These  items  are  discussed  in  appendix  C. 


Distribution 


Copies 


Commandant 

U.S.  Army  Chemical  School 

ATTN:  ATZN-CM-CC  (Mr.  Barnes)  1 

Fort  McClellan,  AL  36205-5020 

NASA  Marshal  Space  Flight  Center 
Deputy  Director 
Space  Science  Laboratory 
Atmospheric  Sciences  Division 

ATTN:  E501  (Dr.  Fichtl)  1 

Huntsville,  AL  35802 

NASA/Marshall  Space  Flight  Center 
Atmospheric  Sciences  Division 

ATTN:  Code  ED-41  1 

Huntsville,  AL  35812 

Deputy  Commander 

U.S.  Army  Strategic  Defense  Command 

ATTN:  CSSD-SL-L  (Dr.  Lilly)  1 

P.O.  Box  1500 
Huntsville,  AL  35807-3801 

Deputy  Commander 

U.S.  Army  Missile  Command 

ATTN:  AMSMI-RD-AC-AD  (Dr.  Peterson)  1 

Redstone  Arsenal,  AL  35898-5242 

Commander 

U.S.  Army  Missile  Command 

ATTN:  AMSMI-RD-DE-SE  (Mr.  Lill,  Jr.)  1 

Redstone  Arsenal,  AL  35898-5245 
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Commander 

U.S.  Army  Missile  Command 

ATTN:  AMSMI-RD-AS-SS  (Mr.  Anderson) 
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